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Abstract 

We study how the propagation speed of interfaces in the Allen-Cahn phase field 
model for phase transformations in solids consisting of the elasticity equations and 
the Allen-Cahn equation depends on two parameters of the model. The two pa¬ 
rameters control the interface energy and the interface width but change also the 
interface speed. To this end we derive an asymptotic expansion of second order for 
the interface speed, called the kinetic relation, and prove that it is uniformly valid in 
both parameters. As a consequence we show that the model error is proportional to 
the interface width divided by the interface energy. We conclude that simulations of 
interfaces with low interface energy based on this model require a very small interface 
width, implying a large numerical effort. Effective simulations thus need adaptive 
mesh refinement or other advanced techniques. 

This version of the paper contains the proofs of Theorem 4.5 and Lemma 5.8, 
which are omitted in the version published in Continuum Mechanics and Thermo¬ 
dynamics. 
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1 Introduction 

In this paper we study how the propagation speed of interfaces in the Allen-Cahn phase 
field model for phase transformations in elastic solids depends on two parameters of the 
model. The model consists of the partial differential equations of linear elasticity coupled 
to the standard Allen-Cahn phase field equation. The two parameters, which we denote 
by fi and A, control the interface energy and the interface width, but variation of these 
parameters also changes the interface speed, or more precisely, the form of the relation, 
which determines the interface speed as a function of the stress field and the curvature of 
the interface. In sharp interface models this relation is called kinetic relation. We use this 
notion also for phase field models. Our goal is therefore to determine the kinetic relation 
for the Allen-Cahn model and the dependence of it on the two model parameters. To this 
end we must derive an asymptotic expansion for the propagation speed of the interface 
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and prove an error estimate for this asymptotic expansion, which holds uniformly in both 
parameters. Our results have consequences for the efficiency of the Allen-Cahn model in 
numerical simulations of interfaces with small interface energy. These consequences are 
also discussed. 

Let n C be a bounded open set with a sufficiently smooth boundary dQ. The 
points of n represent the material points of a solid elastic body. The unknown functions 
in the model are the displacement u{t, x) € M of the material point x at time t, the Cauchy 
stress tensor T{t,x) € 5^, where denotes the set of all symmetric 3 x 3-matrices, and 
the order parameter S{t,x) G M. These unknowns must satisfy the model equations 

-div3,r = b, (1.1) 

T = D{e{V,u)-eS), (1.2) 

dtS = (1.3) 

in the domain [0, oo) x Q. The boundary and initial conditions are 

u{t,x) = U(t, x), it,x) G [0,oo) X dQ, (1.4) 

dng^S{t,x) = 0, (f,x) G [0,oo) X 90, (1.5) 

S'(0, x) = S(x), X G O. (1.6) 

Here b(t, x) G U(t, x) G M^, S(t, x) G M denote given data, the volume force, boundary 
displacement and initial data, dng^ denotes the derivative in direction of the unit normal 
vector nao to the boundary. The deformation gradient S/xu{t, x) is the 3 x 3-matrix of 
hrst order partial derivatives of u with respect to the components x^ of x, and the strain 
tensor 

£{VxU) = -{VxU + (Vxuf) 

is the symmetric part of the deformation gradient, where (Vxu)'^ denotes the transpose 
matrix. The elasticity tensor D : 5^ ^ is a linear symmetric, positive definite 
mapping, e G is a given constant matrix, the transformation strain, and /r > 0 and 
A > 0 are parameters. The elastic energy is given by 

\N{eiVxu),S) = ^(^D{e{Vxu) -eS)) : (e(V,n) -e5), (1.7) 

with the matrix scalar product A : B = Using we obtain for the 

derivative 

9sW(e, S) = -£: D(e(V,,u) - eS) = -e:T. (1.8) 

c > 0 is a given constant and ■0 : M ^ [0, oo) is a double well potential satisfying 

-0(0) = -0(1) = 0, ^(C) > 0 for C / 0,1. 

The precise assumptions on 0, which we need in our investigations, are stated in Theo- 
rem This completes the formulation of the model. 

(HTD and (jl.2l] are the equations of linear elasticity theory. This subsystem is coupled 
to the Allen-Cahn equation (11.31) . which governs the evolution of the order parameter S. 
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The system dLl]) - ([13]) satisfies the second law of thermodynamics. More precisely, the 
Clausius-Duhem inequality is satisfied with the free energy 

r^^{e,S) = \N{e,S) + (1.9) 

From this expression we see that the parameter A determines the energy density of the 
phase interface. We assume that the parameters /U and A vary in intervals (0, /io] and 
(0, Ao], respectively, with /tq > 0 and Aq > 0 chosen sufficiently small. The scaling 
on the right hand side of ()1.3I) is necessary for otherwise the propagation speed of the 
diffuse interface would tend to zero for /r ^ 0 or A —0. 

We give now a slightly sketchy overview of the main results in this article. The precise 
definitions and statements are given in Section [2l and in particular in Section 12.41 

We denote solutions of the Allen-Cahn equation by {uac,T\c, Sac), to distinguish 
them from aproximate solutions, which we construct later. The values SAcit,x) ~ 0 or 
SAcit,x) ^ 1 indicate that at the point x G 0 at time t the crystal structure of the 
material of the solid belongs to phase 1 or to phase 2, respectively. The set of all x € 
with 0 < SAc{t,x) < 1 is the region of the diffuse interface at time t. The level set 

rAc(i) = {a: € I 5'Ac(Ta^) = ^} (1-10) 

belongs to this region. For x € Fac( 1) we denote by SAc{'t,x) G M the normal speed of 
this level set at x, and we call sac the speed of the diffuse interface. For this speed we 
derive an expression of the form 

SAc{t,x) = (soo + A^/^soi) + + A^/^sn) + | In (1.11) 

We call this expression the kinetic relation of the Allen-Cahn model. It is the central 
result of this paper. The remainder term depends on ^ and A, but in Section [6| 
we prove that there exist numbers fiQ > 0 and Aq > 0 and a constant Cs such that the 
L^-norm satisfies 

„<Ce (1.12) 

for all 0 < /r < /io and all 0 < A < Aq, where s^'^^(t) denotes the function x i-)- x) : 

rAc(l) “>■ Therefore (II.lip is an asymptotic expansion for the propagation speed of 
the diffuse interface, which is uniformly valid with respect to the parameters /U and A. 
For sufficiently small the leading term so{t,x) = sqo + A^/^sqi and the second term 
/r^/^si(t, x) = /r^/^(sio+A^/^sii) dominate over the remainder term | In/ip/i 5^"^^ We can 
therefore read off the behavior of the Allen-Cahn model with respect to the parameters 
/i and A from the first two terms in (|l.lll) . 

The terms sqo + and sio + A^/^sn are explicitly given in Theorem 12.31 We 

restrict ourselves here to state the form of the leading term. We assume that 0 < < 

t 2 < 00 are given times. We study the propagation of the interface for t varying in the 
interval [ti,t 2 ]- For t from this interval the leading term is 

so{t, x) = Soo + A^/^soi = ^ - e : (r)(t, x) -h A^/^cifi)r(i, x )^, (1.13) 
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where c is the mobility constant from (jl.3p . ci = y 2i/j{'d) (I'd is computed from the 
double well potential, K,r{t,x) is twice the mean curvature of the surface rAc(i) at the 
point X, and T is the stress field in the solution {t,x) i->- (^u{t,x),T{t,x)) of the trans¬ 
mission problem 


—div3;T 

= b. 

(1.14) 

f 

= D{e{V,u)-eS), 

(1.15) 

N 

= 0, 

(1.16) 

[f]n 

= 0, 

(1.17) 


= m, 

(1.18) 


Here {t,x) S{t,x) is a function, which takes only the values 0 or 1 and jumps across 
the interface 

Tag = {{t,x) E [ti,t 2 ] X H I X E rAc(i)}- 

[lx] and [T] are the jumps of the functions m and T across Tacj and b, U are the volume 
force and boundary data from equations (UdD and (II.4p . The equations (I1.14h . (jl.l5p 
hold in the domain ([ti,t 2 ] x H) \ Tag- For every fixed t the problem (11.141) - (ll.lSp is 
an elliptic transmission problem for the function x !->■ (u{t,x),T{t,x)) in the domain H. 
Therefore t can be considered to be a parameter in this problem. Finally, the expression 

(f) = 

in eauation dl.lSp denotes the mean value of the values T(+) and T^~'> on both sides of 
the interface Tag- 

The terms sio and sn are determined by a more complicated transmission problem, 
for which the transmission conditions are also posed on the interfae Tag; and by a coupled 
system of ordinary differential equations for two functions Sq and S*!, which are needed 
in the construction of an asymptotic solution of the Allen-Cahn model (jl.ip - (jl.51) . The 
coefficients of the second transmission problem and of the system of ordinary differential 
equations depend on the solution (m,T) of (jl.l4p - (I1.18h . Both transmission problems 
together can thus be considered to be a larger transmission problem, which is recursively 
solvable. 

The derivation of the kinetic relation (jl.lip is based on the construction of an asymp¬ 
totic solution for the Allen-Cahn model. From this asymptotic solu¬ 

tion it is seen that the width of the diffuse interface in the Allen-Cahn model is propor¬ 
tional to the parameter 

B = (1.19) 

We call B the interface width parameter. As will be explained in Section [T4l the interface 
energy density is proportional to the parameter 

E = (1.20) 

We call E the interface energy parameter. 

The kinetic relation (II.lip and the equation (I1.19|) together have consequences for the 
efficiency of numerical simulations of interfaces with low interface energy density, which 
we sketch here. A precise discussion is given in Section 12.41 
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The explicit expressions in Theorem 12.81 show that the second term /r^/^(sio + A^/^sii) 
in 111.lip is of a very special form. We therefore argue that this term does not have a 
physical meaning, only the leading term sqo + is physically relevant. This means 

that in (jl.lll) the term 

^(mA) ^ + I 


is a mathematical error term, which in a precise numerical simulation of the evolution 
of the interface mnst be made small by choosing small enongh. Therefore we call 
£’(/^A) model error and F = the error parameter. By (jl.l9p . the interface width 
is proportional to the error parameter F. This means that the interface width depends 
on the size of the model error. The smaller the model error is, which we want to allow, 
the smaller the interface width must be chosen. The total error in a numerical simulation 
consists of the model error and the numerical error. In order to make the numerical error 
small, the grid spacing must be chosen small enough to resolve the transition of the order 
parameter from 0 to 1 across the diffuse interface. When the interface width is small, we 
must therefore choose the grid spacing small, which means that the numerical effort is 
high. 

From (I1.19P and (jl.20p we see that for constant values of the error parameter F the 
interface width is proportional to the interface energy density parameter E. Thus, when 
we want to precisely simulate an interface with small interface energy density, we must 
choose small values for E and E, hence the interface width B = EE becomes very small. 
As a consequence, also the grid spacing must be chosen very small, which means that 
numerical simulations of interfaces with low interface energy based on the Allen-Cahn 
model are not efficient. Of course, the efficiency can be improved by using adaptive mesh 
refinement and other advanced numerical techniques, but still it would be advantageous 
if such tools could be avoided. 

Often the Allen-Cahn model is formulated using the parameters E and F instead of 
H and A. It might therefore be helpful to shortly discuss the form, which our results take 
when this formulation is used. With these parameters the Allen-Cahn equation (jl.8p is 

dtS = -|(5sW(e(V,u),5) + - BA.S)^ 

the free energy (II.9p becomes 


^EBiF S) = W(e, S) + e(^^ + || V,5|2) , 
and the kinetic relation (jl.Iip takes the form 


x) = (^00 + Esqi) + —(<510 + Esii) + 


-'(I) 




From this equation we see that the model error 

B 


B\ 


=—{sio + Esn) + 




is governed by the ratio When one reduces the interface energy density E and one 
wants to keep the model error constant, one must reduce the interface width B by the 
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same proportion. In a simulation of an interface without interface energy, the total model 
error is the sum = Esqi + To reduce the error in such a simulation we must 

reduce E and the fraction hence B must be reduced faster than E. 

The paper is organized as follows. The main results are presented in Section [2l 
where we first state the transmission problems and the system of ordinary differential 
equations, whose solutions are needed to compute the coefficients sqo, • • •, sii in (II.111) . 
These coefficients are explicitly given in Theorem 12.31 Moreover, in this theorem we 
also state properties of the asymptotic solution which is constructed 

in later sections. In particular, we state the scaling law (ll.lQp for the width of the 
diffuse interface. These properties are needed in Section 12.41 where we precisely discuss 
the model error and the numerical efficiency. The estimate (I1.12p . which is the most 
important mathematical result of this paper, is stated in Theorem 12.81 

Sections [3] - [5] contain the proof of Theorem 12.31 In Section [3] we construct the ap¬ 
proximate solution That is, we state the inner and outer expansions 

which define the function In these asymptotic expansions functions 

appear, which are obtained as solutions of systems of algebraic and differential equa¬ 
tions. These systems are also stated in Section O The system for the outer expansion 
can be readily solved, and the solution of the system of ordinary differential equations 
for the inner expansion is more involved and is discussed in Section 01 In two equations 
of this system a linear differential operator appears with kernel different from {0}. In 
order that these differential equations be solvable the right hand sides must satisfy or¬ 
thogonality conditions. The right hand sides contain the coefficients sqo, ■ ■ ■ ,sii of the 
kinetic relation (II.Iip . The orthogonality conditions dictate the form of these coefficients; 
therefrom the equation (jl.Iip results. In Section [5] we verify that ) is 

really an asymptotic solution of the model equations (HI]) - ([L5]) and prove the necessary 
estimates. 

In Section [6] we prove the estimate ()1.12l) . The proof uses the residue, with which the 
function satisfies the equations (jl.ip - (II.3p . The main difficulty in 

the proof is that though we want to prove that is bounded uniformly with respect to 
/r and A, the residue term itself is not bounded for A ^ 0, but instead behaves like 

In the bibliography of [8| we gave many references to the literature on existence, 
uniqueness and asymptotics for models containing the Allen-Cahn and Cahn-Hilliard 
equations. We refer the reader to that bibliography and discuss here only some publica¬ 
tions, which are of interest in the construction of asymptotic solutions. 

We believe that for the model (II.ip - ()I.3|) an asymptotic solution was constructed 
and used to identify the associated sharp interface problem for the first time in |23j . 
following earlier such investigations for other phase field models. For example, in m 
these investigations were carried out for a model from solidification theory, which consists 
ot the Allen-Cahn equation coupled to the heat equation. 

The considerations in mm are formal, since it is not shown that the asymptotic 
solution converges to an exact solution of the model equations for /r ^ 0. Under the 
assumption that the associated sharp interface problems have smooth solutions, this was 
proved in [28| for the Allen-Cahn equation, in [TO] for the Cahn-Hilliard equation, in m 
for the model from solidification theory and in [1] for a model consisting of the Cahn- 
Hilliard equation coupled with the elasticity equations. The proofs use variants of a 
spectral estimate derived in |16] . For the model from solidification theory the associated 
sharp interface model is the Mullins-Sekerka model with surface tension. 


6 
















In [15] an asymptotic solution for the Cahn-Hilliard equation has recently been con¬ 
structed with a method different from the one used in ini, and which is similar to our 
method. 

In |25j the numerical efficiency of simulations based on the phase field model consisting 
of the Allen-Cahn equation coupled to the heat equation is studied. It is shown that for 
suitably chosen coefficients of the model the second order term in an asymptotic expansion 
of the solution vanishes. By arguments similar to the ones we gave in the above discussion 
it is seen that this improves the numerical efficiency of the model. This result has been 
improved and generalized in PEaEU. A similar idea is also present in m- 

Since the construction of asymptotic solutions is based on sharp interface problems, 
a rigorous analysis of these problems is of special interest. Of particular interest is 
the Hele-Shaw problem with surface tension, since this is the sharp interface problem 
associated with the Cahn-Hilliard equation. Existence, uniqueness, and regularity of 
classical solutions of this problem have been investigated in [181 119] 120) . In |2T| it is 
shown that if the initial data are close to a sphere then a classical solution exists and 
converges to spheres. Existence of solutions to the Mullins-Sekerka problem mentioned 
above has been shown in m- 

The model (jl.ll) - (II. 6p describes the evolution of phase transitions in a solid when 
temperature effects are negligible. This model is the prototype of a large class of models 
obtained by extensions and generalizations of the model, which are used in the engineering 
sciences to simulate the behavior of complex and functional materials. From the very 
large literature in this field we cite here only [IIl|29l[3ni|3ll|32]. 

2 The kinetic relation 
2.1 Notations 

For given fixed times 0 < < ^2 < oo let 

Q = X H C 

The construction of the asymptotic solution ) is based on a surface 

r(A‘^)(i), which for ti < t < t2 moves in H and which will be the level set 

= {x G H I = ^}. 

We set 

= {{t,x) e Q \ X e (t)}. (2.1) 

To simplify the notation we often drop the index A or both indices n and A and write 
and T^^) or simply r(t) and T. Similarly, we often write for the 

asymptotic solution and use the same convention also in other notations. Both indices 
are specified if the dependence on A becomes important. 

The precise definition of the family is given in the nect section, 

and in Section [23] we associate with the level set = rAc(i) introduced 

in (jl.IOp . To introduce notations we assume here that T is a known, orientable, three 
dimensional C'*^-manifold with k > 1 sufficiently large embedded in Q such that T{t) is 
a regular two dimensional surface in H for every t G [ti,t 2 ]- Let 

n : r ^ (2.2) 
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be a continuous vector field such that n{t,x) € is a unit normal vector to r(t) at 
X G r(t), for every t € [ti,t 2 ]. For 5 > 0 and t € [ti,t 2 ] define the sets 

Us{t) = {x G n I dist(x,r(t)) < (5} and Us = {{t,x) G Q | x ^Us{t)}. (2.3) 

We assume that there is 5 > 0 such that Us C Q. Since F is a regular C^-manifold in Q, 
then 5 can be chosen sufficiently small such that for all t G [^ 1 ,^ 2 ] the mapping 

{rj, 0 '-t x{t, ri,^) = rj + ^n{t, t]) : F(t) x (-6, <5) ^ Us{t) (2.4) 


is bijective. We say that this mapping defines new coordinates {r],^) in Us{t) and (t, rj, ^) 
in Us- If no confusion is possible we switch freely between the coordinates {t,x) and 
In particular, if (t,x) i->- w{t,x) is a function defined on Us we write 
for w[t,x{t,r],^)'j, as usual. 

We use the standard convention and denote for a function w defined on a subset U 
of Q by w{t) the function x i-)> w(t,x), which is defined on the set {x | (t,x) G U} C 
If rc is a function dehned on Us{t) \T(t), we set for r/ G F(t) 


(rj) 

idnwY~\v) 

[w]{rj) 

[dnw]{r]) 

{w){r]) 


lim w(^r] ± |n(t, r/)), 
e>o 

(9* 

lim ^u;(r/ + ^n(t,r/)), f G N, 
e>o ^ 


Ql 

hm TT-rtc 

dU 

e<o ^ 


{r] + Cn{t,rj)), 


w^^\r]) — \r]), 


i G N, 


(5»W(r/) - ((9»( \ri), 

^{w^^^Tj) + w^~\rj)), 


provided that the one-sided limits in these equations exist. If w is defined on \ F, we 
set 

w^^\t,r]) = {w^^\t)){rj), (0»(=^)(t,r/) = {rj), 

and define [w]{t,r]), {w){t,r]), [dn'w]{t,r]) as above. Let Ti(r/), T 2 {r]) G be two orthog¬ 
onal unit vectors to T(t) at r] G F(t). For functions w : T{t) ^ W, W : F(t) ^ we 
define the surface gradients by 


Vrre = {driw)Ti + {dr2w)T2, (2.5) 

VrIF = {dr,W)0Tl + {dr,W)(^T2, (2.6) 

where for vectors c, d G a 3 x 3-matrix is defined by 

C®d= (Cjdj)ij=p2,3- 

With (j2.5h . (12.6h we have for functions w :Us{t) and IF : Us{t) —)> at t/ G T{t) 

the decompositions 


VxW = {dnw)n + Vvw, 
V^IF = (9nIF) 0 n + VrIF, 
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(2.7) 

( 2 . 8 ) 



where n = n{t,r]) is the unit normal vector to T(t). 

The normal speed of the family of surfaces t i—>• r(t) is of fundamental importance in 
this paper. Therefore we give a precise definition. 


Definition 2.1 Let m{t,ri) = (^m'{t,r]),m"{t,r])'j € M x be a normal vector to T at 
{t, rj) E r. The normal speed of the family of surfaces 1 1-7> T{t) at rj ^ r(t) is defined by 


s{t,r]) 


—m'{t, r]) 
m"{t, rj) ■ n{t, rj) ’ 


(2.9) 


with the unit normal vector n{t,r]) E to T(t). 


Note that with this definition the speed is measured positive in the direction of the normal 
vector field n. Since m"{t,rj) E is a normal vector to r(t), the denominator in ()2.9p 
is different from zero. 

\i uj = [oj',oj") E M X is a tangential vector to T at (t,?/) with oj' 0, then with 

the unit normal n(t,r]) E to T{t) the vector {—u" ■ n,uj'n) is a normal vector to T at 

hence ()2.9I) implies that the normal speed at r/ E r(t) is given by 


s{t,ri) 


n ■ UJ 


n 


uj'n ■ n 


n ■ UJ 


n 


uj' 


For later use we prove the following 


( 2 . 10 ) 


Lemma 2.2 Let x E ^ 5 (^ 0 ) a point having the representation x = rj + n{t,r])f in 
the {rj,fj)-coordinates, where ij = r]{t,x) E r(t) and ^ = f,(t,x). Then the normal speed 
satisfies 

s{to,r]) = n{to,ri) ■ dtr]{to,x) = -dtf{to,x). ( 2 . 11 ) 

The tangential component of the vector dtr]{tQ,x) E to the surface r(to) is equal to 
-^dtn{to,r]{to,x)). 

Proof: By definition of Us, there is a neighborhood U of to in [^ 11 ^ 2 ] such that {x} xU C 
Us, which implies that x has the representation 

X = T/(t, x) + f{t, x)n{t, rj{t, x)). 

for all t E t/. We differentiate this equation and obtain 

0 = dtx = ndt^ + ^ dtn + dtp. (2-12) 

From 0 = dfl = i9t|np = 2n • dtn we see that dtn is tangential to r(t), hence (j2.12ll 
implies that the tangential component of dtp is equal to —^dtn. Multiplication of (j2.12p 
with n yields 

dtf = -n-dtp. (2.13) 

Since dt{t, p{t, x)) = {l,dtp{t, x)) is a tangential vector to F, it follows from (I2.10p that 
s = which together with (j2.13p implies (j2.1ip . ■ 
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2.2 The evolution problem for the level set 

The level set T = of defined in () 2 . 1 [) is determined by an evolution problem for 
the family of surfaces t i—)> r(t). To state this evolution problem let Af be the operator, 
which assigns the normal speed to the family 1 1 —>■ r(t), i.e. 

s{t, x) = Af{r){t, x), 

with s{t,x) = s^^\t,x) defined by (12.91) . The evolution problem is given by 

AA(r)(t) = /c(^)(r(t)), (2.14) 

where is the non-local evolution operator, which has the form 

(r(t)) (x) = so (r, Kr, Ai/2) (t, x) + {u, f, T, So, Si, A^/^) (t, x), (2.15) 

for X € r(t). Here {u,T,u,T, Sq, Si) is the solution of a transmission-boundary value 
problem for a coupled system of elliptic partial differential equations and ordinary dif¬ 
ferential equations, which can be solved recursively. Krit,x) denotes twice the mean 
curvature of the surface r(t) at x € r(t). With the principle curvatures ki, K 2 of r(f) at 
X G r(t) we thus have 

Kr{t,x) = Ki(t, x) -I- K2it,x). 

The transmission condition is posed on r(t). Therefore the functions u, T, u, T and Si 
depend on r(t). We first state and discuss the transmission-boundary value problem. 
The precise form of the functions sq •Si is given in Theorem 12.81 following below. 

Let S'. Q \ r —>• {0,1} be a piecewise constant function, which only takes the values 
0 and 1 with a jump across L. The sets 

7 = {(t, x) € Q \ r I S(t, x) = 0 }, 7(t) = {x € H \ r(f) I (f, x) G 7}, 

7 ' = {{t,x) e Q\T \ s{t,x) = 1}, j'{t) = {x G H\ r{t) I (t,x) G 7 '} 

yield partitions Q = 7 U L U 7 ' and H = 7 (f) U r(f) U 7 '(t) of Q and H, respectively. If 

X belongs to 7 (f) or 7 ^(t), then the crystal structure at the material point x at time t 

belongs to phase 1 or phase 2, respectively. We assume that the normal vector field n 
given in (|2.2p is such that the vector n{t,x) points into the set ^'{t) for every x G r(t). 

The transmission-boundary value problem can be separated into two transmission¬ 
boundary value problems for the elasticity equations and a boundary value problem for 
a coupled system of two ordinary differential equations. To state the complete problem 
we fix t G and assume that r(t) is known. In the first transmission-boundary 

problem the unknowns are the displacement x 1 —)• {t(t, x) G and the stress tensor 
x ^ T{t,x) G <S^, which must satisfy the equations 


-div3;T 

= b. 

(2.16) 

f 

= D{e(y,u)-eS), 

(2.17) 

N 

= 0, 

(2.18) 

[f]n 

= 0, 

(2.19) 


= U(f), 

(2.20) 
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with b and U given in (jl.lj) and (11.41) . In the second transmission-boundary problem the 
unknowns are the displacement x i-)- u{t, x) € and the stress tensor x !->■ T{t, x) € 
and the problem is 


-diVa,T 

= 0, 

T :e \ 

(2.21) 

t 

= D(^e{Vxu) 

- e -—^ 

V’"(5)^ 

(2.22) 

[u] 

= 0, 


(2.23) 

[f]n 

= 0, 


(2.24) 

"(‘)U 

= 0. 


(2.25) 


The equations (I2.16D . ()2.17p and (I2.2ip . (12.221) must hold on the set n\r(t), whereas the 
equations ()2.18l) . p2.19l) and (12.2311 . (I2.24p are posed on r(t). 

In the bonndary valne problem for the ordinary differential equations the unknowns 
are S'o : M —>■ M, Si : T x M —>■ M and sq : T —>■ M. We use the notations r/, C) = 
d(^Si{t,r]X), Si{t,r]X) = In this problem not only t, but also rj G r(t) is a 

parameter. For all C € IR and all values of the parameter r] G r(t) the unknowns must 
satisfy the coupled ordinary differential equations 

^'(So(C))-5^'(C) = 0, (2.26) 

^''{So{0)Si{t,ri,C)-SUt,vX) = (2.27) 


and the boundary conditions 


'S'o(O) = lim So(C) 

Z (—>—00 

lim Si{t,r],C) 
(—>—00 

lim Si{t,r],C) 

(->+oo 

Si{t,r],0) 


0, lim So(C) = 1, 

(—>oo 

^"( 0 ) ’ 

r{l) ’ 

0 , 


with the right hand side of p2.27p given by 


Fi (t, v,C)=e: ([f ] it, r?)So(C) + (t, v)) 


+ 


soit^v) 


xF^Krit,'n)'^S'oiC), 


(2.28) 

(2.29) 

(2.30) 

(2.31) 


(2.32) 


where the constant c > 0 is given in (II.3p . 

The linear elliptic system (12.1611 , p2.17l) differs from the standard elasticity system only 
by the term —DeS. This term is known since r(t) is given. Under suitable regularity 
assumptions for the given fnnctions b and U and very mild assumptions on the regnlarity 
of the interface r(t) the problem has a unique weak solution [u,T). This can be proved 
by standard methods from functional analysis. Of course, the regularity of the solution 
depends on the regularity of b, U and r(t). 
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After insertion of the stress tensor T from this solution into (j2.22p . the equations 
(j2.21l) - (j2.25l) form a transmission-boundary value problem of the same type as (12.161) - 
(j2.2np . with unique solution {u,T) determined by the same methods. 

We also insert T into (j2.29p . (I2.30p and (|2.32l) . which determines the right hand 
side of the differential equation ()2.27p and the boundary conditions (j2.29p . (12.301) posed 
at ±oo. The nonlinear differential equation (I2.26P has a unique solution Sq satisfying 
the boundary conditions (I2.28h . By insertion of Sq into (I2.27h and (j2.32p . equation 
(j2.27p becomes a linear differential equation for Si, however with an additional unknown 
function sq in the right hand side. This function is constant with respect to C- We sketch 
here the procedure used to determine sq- This procedure is standard in investigations of 
the asymptotics of phase field models: 

The second order differential operator is selfadjoint in the Hilbert 

space L^(M) with a one dimensional kernel spanned by the function Sq. This is seen 
by differentiating the equation (I2.26p . From functional analysis we thus know that for 
Fi € L^(M) the differential equation ('^"(S'o) — d‘^)w = Fi has a solution w G L^(R) if 
and only if the orthogonality condition 

/ OO 

Fiit,v,C)S'oiOdC = 0 (2.33) 

-OO 

holds. It turns out that though the function Fi defined in (I2.32|) does not in general 
belong to L^(]R) and the solution Si{t,r],-) is not sought in which is seen from 

the boundary conditions ()2.29p , (j2.30p , the orthogonality condition (I2.33h is sufficient for 
the solution to exist. Comparison with (I2.32p shows that (|2.33l) can be satisfied by 
choosing the constant so{t,r]) suitably. This defines the function sq : T ^ R uniquely. 
Since Fi depends on T, Kr and it follows that also sq is a function of these variables: 

so{t, rj) = so{f, Kr, {t, x). 

The explicit expression for sq obtained in this way is stated below in ()2.40p . In fact, 
so(T', Kr, A^/^) is the first term on the right hand side in the expression (j2.I5p for 

(r(t)). 

The procedure sketched here is discussed precisely in Section 14.21 when we determine 
the second term si in ()2.15p . which is obtained from a similar, but more complicated 
boundary value problem. 

2.3 The asymptotic solution and the kinetic relation 

To state the properties of the asymptotic solution and the kinetic relation in Theorem l2.3l 
we introduce some definitions. 

We need in our investigations that the second derivatives and of the 

double well potential at the minima 0 and 1 are positive, and we set 

a = min I (0), ?^"(I)|. 

Depending on the parameters A and /x, we partition Q into the inner neighborhood 
of T, into the matching region Q^^tch outer region Qout^. These sets are 
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defined by 


^(mA) r,, 3(^A)i/2|ln/i|-, 

Qinn = 1^1 < 2 - - -/’ 

oLl’ch = {(*.».«) 


3(MA)V^j^ 300)^1^1 

2 a - 1^1 - a 


/’ 


(2.34) 




We always assume that the parameters A and /r satisfy 0 < A < Aq and 0 < fi < fiQ, 
where Aq, fJ-o are fixed constants satisfying 


^ _2 3(/ioAo)^/^|ln/io| 

Ho<e , - < d. 

a 

The first condition is imposed for purely technical reasons and guarantees that the func¬ 
tion n !->■ In^ul is increasing, the second condition guarantees that Qmatch ^<5 

and that n is a nonempty, relatively open subset of 14^. 

By (I2.18j) . the function ft : Q —is continuous at every point (t, ??) € T, but the 
first and higher derivatives of u in the direction of the normal vector n(t, x) can jump 
across T. For these jumps we write 


We also set 


u*{t,r]) = [dnu]{t,r]), 

(2.35) 

a*{t,rj) = [dlu\{t,r]). 

(2.36) 

Cl = f \J2'ip{'d) d'd. 

Jo 

(2.37) 


Theorem 2.3 Suppose that the double well potential ijj G C'^(M) satisfies 


fi{r) > 0, for 0 < r < 1, 

iP[r) =fi'{r) =0, for r = 0,1, (2.38) 

a = min | '^"(O), '0"(1)| > 0. 

Moreover, suppose that satisfies the symmetry condition 

^(^-C) = ^(^ + C), CgK. (2.39) 

Assume that there is a solution T of the evolution problem (j2.14p . ()2.15p with sq = 
so(r,Kr,AV2) :r^]R given by 

so = ^(-e: (r) + Ai/2^iKr), (2.40) 

and with si = si('u, T, T, So, Si, A^/^) : T —)> M defined by 

Si = Sio + Ai/2g^^ = ^ So, Si) + So), (2.41) 
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where 


•sio 


•Sll 


1 1*00 -I POO 

+ -e : (f) / dC + - r'{So)SfS'o dQ 

J—oo ^ J—oo 

c 

- e:De{a*^n + Vru*) So{C)So{-C) dC- 

J—oo 


(2.42) 

(2.43) 


In (I2.40p and (12.421) . (I2.43P we have Sq = 5'o(C) and Si = Si{t,r]X), for all other 
functions the argument is ft, rj). The positive constant c is defined in (|1.3p . The notations 
[•] and {■) are introduced in Section ro In particular, we have 


/ e:f \ Ir e: f(+) e 
\V>"(5)/ 2 V ^"(1) ^"(0)^' 

With these funetions the normal speed s{t,r]) ofT{t) at rj £ T{t) is thus given by 

s{t,r]) = so{t,rj) + ji^/‘^si{t,rj,X^/‘^) = so{t,7]) + p^/‘^{sio{t,rj) + X^/'^snft,!])). (2.44) 

We assume moreover that the solution T is a -manifold and that the functions u and 
u defined by the evolution problem satisfy u G C^('y U u E ( 7^(7 U 7 ', M^) and that 

u has -extensions, u has -extensions from 7 to 7 U F and from 7 ' to 7 ' U F. For 
the given right hand side of o we assume that b E C^{Q). 

Under these assumptions there is an approximate solution {u^^\T^^\ of the 

Allen-Cahn model dn]) - dLSp, for which F is the level set 


F = |(t,x) E <3 1 


(2.45) 

and which satisfies the equations 



-div^r(^) = 

= b + fi^^\ 

(2.46) 

r(^) = Zl(e(V, 


(2.47) 

{jlX)2 ^ 

P'2 ' 

(2.48) 

II 

cz 

(t, x) E [tl,t2] X dl}, 

(2.49) 


(t,x) E [tl,t2] X dl}, 

(2.50) 

where to the right hand sides 

there exist nonnegative constants Ki, 

...,iF5 

such that for all p E (0, po] and all X E (0, Aq] 


II II 

ll/i 

- (0 "^1 ’ 

(2.51) 

II II 

lUi 

< U^K2, 

(2.52) 

II II 

11/2 

< |ln/x|2(0"iF3, 

(2.53) 

II II 

IIJ 2 


(2.54) 

II f II 

II /3 IIl°°(90) 

< 

(2.55) 
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In the neighborhood ^ of V the order parameter in the approximate solution is of the 
form 

where the monotonieally inereasing transition profile So : M —)• M and the function Si : 
r X M —)• M are given as solution of the eoupled problem (12.261) - (12.321) . and where 
S2 : r X M ^ M satisfyies 82 ( 1 , r/, 0) = 0 and 

|S2(t,??,C)l < (^(1 + ICI), /or (i, T/, C) e r X M, (2.57) 

with a constant C independent of {t,r]X)- 

We mention that the positive constant c in ()1.3[) does not play a major role in the analysis 
and could be replaced by 1. We refrain from replacing it to show how c appears in the 
kinetic relation. 

The proof of this theorem forms the content of Sections El-El We remark that 
the symmetry assumption (j2.39p for the double well potential if serves to simplify the 
computations in the derivation of the asymptotic solution. Without this assumption the 
term si in the kinetic relation (|2.44p would contain other terms in addition to the terms 
sio and sn given in (j2.42p an (j2.43p . 

The regularity properties of T and of u, u are of course not independent, since u 
and u are solutions of the elliptic transmission problems (I2.16|) - (|2.20p and (|2.2ip - 
(|2.25l) . respectively. Therefore the regularity theory of elliptic equations shows that u 
and u automatically have the differentiability properties assumed in the theorem if the 
manifold T and the right hand side b are sufficiently smooth. 

Since by definition of and (I2.34h we have 

meas(g[^^^ U ^\, 

we immediately obtain from ()2.51l) - (j2.54p the following 

Corollary 2.4 There are eonstants Kq, Ky such that for all 0 < p < fiQ and all 0 < 

A < Ao 

ll/r’lliHO) S ItafiVA'e, (2.58) 

|/7"’llt‘«3) S (2-59) 

The leading term sq given in (I2.40p can be written in a more common and more general 
form. To give this form, we need a result on the jump of the Eshelby tensor. The Eshelby 
tensor to the solution (tt,T) of the transmission problem (I2.16p - (j2.20p is defined by 

C(V,u, S) = V’/.(e(V,u), S)I-{I + V^uff, (2.60) 

where I G 5^ is the unit matrix and where 

V>4e,5)=W(6,5) + ^^(5) (2.61) 

is that part of the free energy if"^ defined in (II.9|) without gradient term. The last term 
on the right hand side of (|2.60p is a matrix product. We use the standard convention to 
denote the matrix product of two matrices A G and B G by AB G 
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Lemma 2.5 Let {u,T) be the solution of the transmission problem (j2.16p - (j2.20p and 
let n be a unit normal vector field to T{t). Then the jump [C] of the Eshelby tensor to 
{u, T) across F satisfies 

n-[C]n = ^mS)]-e:{f). (2.62) 

This result is known [2]. In [7] it is stated as equation (3.4) and proved on pages 154, 
155. 

Corollary 2.6 The leading term sq of the kinetic relation defined in ()2.40p satisfies 

So = — (n ■ [C]n + X^^'^ciKr^. (2.63) 

This corollary follows immediately from ()2.62p . since by assumption ()2.38l) we have 
[?^(5)] = 'i/>(l) — ip{0) = 0, which implies that n • [C]n = —e ; (T). ■ 

2.4 Consequences for numerical simulations 

In this section we discuss the consequences of Theorem 12.31 for numerical simulations of 
interfaces with small interface energy. 

In many functional materials the phase interfaces consist only of a few atomic lay¬ 
ers. For interfaces with such small width mathematical models with sharp interface are 
appropriate. We therefore base the following considerations on the hypothesis that the 
propagation speed of the interface in the sharp interface model is a good approximation 
to the propagation speed of the interface in the real material. The model error of the 
Allen-Cahn model is then the difference of the propagation speed of the sharp interface 
and the propagation speed of the diffuse interface in the phase field model. The param¬ 
eters pL and A in the Allen-Cahn model should be chosen such that this model error is 
small and such that numerical simulations based on the Allen-Cahn model are effective. 

To make this precise we must first determine the sharp interface model to be used. 
The model consists of the transmission problem (|2.16p - (|2.20l) combined with a kinetic 
relation. To find this relation, one proceeds in the usual way and uses that by the second 
law of thermodynamics the Clausius-Duhem inequality 

dt'f’sharp T diVj; ^sharp ^ ‘ b 

must be satisfied to impose restrictions on the form of the kinetic relation. Here V’sharp 
denotes the free energy in the sharp interface problem and ^sharp is the flux of the free 
energy. We use the standard free energy and flux 

Aharp{e{V,u), S) = W(e(V,u), S) + X^/^ci [ da, (2.64) 

Jrit) 

qsharp{T,S) = - f -Ut. 

The last term on the right hand side of (|2.64l) is the interface energy, hence is 

the interface energy density. It is well known that if ('u(t),T(t)) is a solution of the 
transmission problem (I2.16P - ()2.20p at time t and if the interface T{t) in this problem 
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moves with the given normal speed Ssharp(i, a^) at x G r(t), then the Clausius-Duhem 
inequality holds if and only if the inequality 


■Ssharp(i, x) (^n{t, x) ■ [C]{t, x)n{t, x) + A^/^ciKr(i, x)^ > 0 (2.65) 

is satisfied at every point x G r(t). A proof of this well known result is given in [3], 
however only for the case where A = 0 in (j2.64p . The proof can be readily generalized to 
the case A > 0. 

A simple linear kinetic relation, for which (I2.65h obviously holds, is 


^sharp 


— (n • [Cln + A^/^ciKr) • 

Cl 


( 2 . 66 ) 


The sharp interface problem thus consists of the transmission problem (I2.16P - (12.201) 
combined with the kinetic relation (j2.66p . For this problem the Clausius-Duhem inequal¬ 
ity is satisfied. 

We can now define the model error. To this end let 5^^^) be the asymp¬ 

totic solution in the domain Q = [ti,t 2 ] x Q constructed in Theorem 12.31 where by 
(I2.45p . the manifold T is the level set = i}. Let i G [^ 1 ,^ 2 ] be a hxed time and 

let(u^(^, 5^) be the exact solution of the Allen-Cahn model (|l.ip - (jl.3p in the 

domain [t,t 2 ] x D, which satishes the boundary and initial conditions 


U'^l{t,x) = U(Lx), 

{t, x) G [i, t 2 \ X dCt, 

(2.67) 

^ngaSAci^’X) = f!^^^\t,x), {t,x) G [i,t 2 ] X dn, 

( 2 . 68 ) 

Ci& 

II 

’fk 

), X G D, 

(2.69) 

is the rieht hand side of (I2.5UD. 

The level set of the order parameter 5^ is 

Tag = |(t,x) g Q 


(2.70) 


denoted by 


Let Tsharp C Q be the sharp interface in the solution of the sharp interface problem (j2.16p 
- (j 2 . 20 p . (12.661) . which satisfies the initial condition 


^sha,rp{i) = T{i). 


(2.71) 


The normal speeds of the different surfaces are 

S = = W(r('^^)), SAC = ^^sharp = , 

where M is the normal speed operator introduced at the beginning of Section 12.21 Of 
course, Ssharp is given by (j2.66p . Note that the functions SacH^), Ssharp(i) are 

defined on the same set, since the initial condition (|2.69l) and (|2.7ip together imply 
^Ac{i) = r(f) = Tsharp (t)- 

Definition 2.7 We call the function £ = : r(t) ^ M defined by 

£ = SAc(i) - Ssharp(i) (2.72) 

the model error of the Allen-Cahn model at time t to the parameters p, and A. 
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We can now discuss the choice of the parameters /x and A. Since (I2.66|) coincides with 
the leading term sq in the asymptotic expansion (|2.44l) of the kinetic relation of the 
Allen-Cahn model, which is seen from (j2.63p . we have 


■Ssharp — 'SQ- (2.73) 

Therefore (12.441) yields 

£ = 'SAC“'Ssharp = 'SAC“'So = (sAC “•s) + (s — Sq) = (sAC “'S) + /.i^^^(sio + A^'^^Sll). (2.74) 

The difference sac “ s between the propagation speeds of the exact solution and the 
asymptotic solution tends to zero for /x ^ 0 faster than the term /x^/^sio, the 
convergence is uniform with respect to A. This is the basic result, which allows to discuss 
the optimal choice of /x and A. The precise result is 

Theorem 2.8 There is a constant Cg > 0 such that for all 0 < fi < fiQ and all 0 < X < 
Aq we have the estimate 


||sAc(t) - 'S(t)||i2(r(f)) < Qlln/xp/x- (2.75) 

The proof of this theorem is given in Section O 
(j2.74p and (j2.75p together yield 

(2.76) 

with a constant C, which can be chosen independently of A. By this inequality, /x^/^ 
controls the model error. Therefore we write F = /x^/^ and call F the error parameter. 
Moreover, since A^/^ci is the interface energy density, we call E = A^/^ the interface 
energy parameter. Also, since by (I2.56P the interface width is proportional to (/uA)^/^, 
we call B = (^A)^/^ the interface width parameter. These three parameters and the 
propagation speed sac are connected by the fundamental relations 

B = EF, (2.77) 

SAC = — n ■ [C]n + ckyE + £[E, F], (2.78) 

Cl 

ll^[^,^]|lL2(r(,'))<CF, (2.79) 

where we use the notation S[E,F] = The first equation is an immediate conse¬ 

quence of the definition of the parameters, the second is obtained by insertion of (I2.66p 
into (j2.72p . and the last inequality is just a restatement of (I2.76p . 

Now assume that we want to use a phase field model to numerically simulate the 
propagation of a phase interface. In such a simulation the numerical effort is proportional 
to h~P, where h denotes the grid spacing and where the power p > 1 depends on whether 
we want to simulate a problem in 2-d or in 3-d and it depends on the numerical scheme 
we use. In order for the simulation to be precise, we must guarantee that the model 
error and the numerical error are small. To make the numerical error small, we must 
choose the grid spacing h small enough to resolve the transition of the order parameter 
across the interface, which means that we must choose h < B, hence we have h~^ > B~p. 
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Therefore we see that the numerical effort of a simulation based on a phase field model 
is measured by the number B~p. We call the number 


the parameter of numerical effort. For a simulation based on the Allen-Cahn model we 
see from (j2.77p that the numerical effort is 

enum = {EF)-P. 

Assume that the interface, which we want to simulate with the Allen-Cahn model, has 
very small interface energy density. Such interfaces are common in metallic or functional 
materials. For such materials the interface energy parameter E is small. To make the 
model error small, we must also choose the error parameter F small, which means that 
the numerical effort parameter Cnum = {EF)~'p is very large as a product of two large 
numbers E~P and F~'p. 

To be more specific, we consider an interface without interface energy, which means 
that the free energy V'sharp does not contain the last term on the right hand side of (|2.64p . 
From (|2.66p we see that the propagation speed of the sharp interface with zero interface 
energy density is 

Ssharp ~ ^ ■ [CJri. 

Cl 

From this equation and from (I2.78P we see that in this case the total model error, which 
we denote by £^totai, is 


^total Sac Ssharp CKy'E + £[E, F]. 

This means that the term cKrE is now part of the total model error. This term does not 
vanish identically, since we cannot set A = 0 in the Allen-Cahn equation (jl.3p . Instead 
the values of A and of F = A^/^ must be positive. 

If we prescribe the F^-norm £ 1^2 = ||^^totai|lL 2 (r(f)) of fh® total model error, we must 
therefore choose the parameters E and F such that 

c||Kr||L2(r(f))F-I-||£’[F,F]||^2(r(f)) < £l^, (2.80) 

EE = max, (2-81) 

where the second condition is imposed by the requirement to make the numerical effort 
Cnum = {EF)~P as small as possible. To discuss this optimization problem, we assume 
first that the term sio in the asymptotic expansion (I2.44p of the kinetic relation of the 
Allen-Cahn model is not identically equal to zero. In this case we conclude from (I2.74p 
and (j2.75p by the inverse triangle inequality that for sufficiently small A^/^ = E and for 
sufficiently small = F 

+ ('SAC - 'S)|lL2(r(£)) 

> Ai^^^l|sio|lL2(r(t)) “ (F-^)^'^^l|sii||i;,2(r(£)) - ||sac - 'S|lL2(r(t)) 

> F^^^(l|sio|lL2(r(£)) - A^/2||sii||^2(r(f)) -Q|ln/r|V^^^) 

> F^'^^(ll'Sio|lL2(r(£)) - 2 ll'®iollL2(r(f))) = 2 
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This inequality and (12.801) imply that the solution {E, F) of the optimization problem 
(j2.80p . (12.811) satisfies 


F < 


l'®iollL2(r(f)) 

From this result we obtain 




l'5io|lL2(r(t)) 


Sl2 


and 


E < 


c||'«r|lL2(r(£)) 


£ 12 . 


Corollary 2.9 Let Tmax denote the total model error of the Allen-Cahn model in the 
simulation of an interface without interfaee energy. If the term sio in the asymptotie 
expansion (I2.44h of the kinetic relation of the Allen-Cahn model is not identieally equal 
to zero, then the interface width B satisfies 


B = EF < 


2 

c||'Sio|lL2(r(t))ll^r|li2('r(£)) 



(2.82) 


In a numerical simulation of an interface without interface energy the parameter of nu¬ 
merical effort satisfies 


Cnum ^ 


'5iollL2(rd))ll^r|lL2(r(£)) Y 


2 T|2 




(2.83) 


with a power p > 1 depending on the space dimension and the numerieal method used. 


The interface width thus decreases with the square of the model error. Since the time step 
in a simulation must be decreased when the grid spacing h in x-direction is decreased, 
the number p can be larger than 4 in a three dimensional simulation. From (12.831) we thus 
see that the numerical effort grows very rapidely when the required accuracy is increased. 
The Allen-Cahn model is therefore ineffective when used to accurately simulate interfaces 
with low interface energy. 

If the term sio vanishes identically, then the same considerations show that instead of 

(I2.82P and (I2.83P we would have B = 0{£^2^) and enum > C£^ 2 ^. The numerical effort 
would still grow fast when the required accuracy is increased, though less fast than for 
sio 7 ^ 0. However, a close investigation of the terms in the definition (j2.42p of sio, which 
we do not present here, shows that only in very exceptional situations one can expect 
that sio vanishes identically. 

In Corollary 12.91 we assumed that the mesh is globally refined. Of course, one can 
improve the effectivity of simulations by using local mesh refinement in the neighbor¬ 
hood of the interface. We do not discuss this question of numerical analysis here, but 
Corollary 12.91 in fact shows that adaptive mesh refinement and other advanced numerical 
techniques are needed to make precise simulations of interfaces with small energy based 
on the Allen-Cahn model effective. 


Comparison to the hybrid phase field model With Corollary 12.91 we can refine the 
comparison given in [ 8 ] of the Allen-Cahn model and another phase field model, which we 
call the hybrid model. The hybrid model was introduced and discussed in [HllllSlElIH]. 
By formal construction of asymptotic solutions we showed in [ 8 ] the following result: 

Let Bp,c(£]^ 2 ) and i?hyb(^L 2 ) be the interface widths in the Allen-Cahn model and the 
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hybrid model, respectively, which result when the model parameters are adjusted to 
model an interface without interface energy with the total model error £ 12 . Then we 
have for £’^2 —>■ 0 that 


Bhyh{£L^) = 0{£l2), Bac{£l^) = o{1)0{£l2) = o{£l2). 


The Landau symbol o(l) denotes terms, which tend to zero for Tmav —)• 0. This result 
was obtained under the assumption that estimates corresponding to (j2.75l) hold true for 
both models, without proving this assumption. 

To achieve a prescribed small value £ 1^2 of the total model error we must therefore 
choose the interface width in the Allen-Cahn model smaller than in the hybrid model. 
Consequently, the hybrid model is numerically more effective, but how much more de¬ 
pends on the rate of decay of the o{£]^ 2 ) term in the result for the Allen-Cahn model. In 
[8] we could not determine this decay rate, since the asymptotic solution constructed in 
[8] for the Allen-Cahn model was only of first order. 

Corollary 12.91 yields this decay rate. Prom the result for the hybrid model and from 
Corollary 12.91 we thus obtain for the parameters Cnum and of the hybrid model and 
the Allen-Cahn model, respectively, that 


„hyb 


< r<f~P 

'^num — ’ 


.AC 


®num — £'£ l ^ 


-2p 


which shows that when the prescribed error £^2 is small, the hybrid model can be quite 
considerably more effectice in numerical simulations of interfaces with low interface energy 
or no interface energy than the Allen-Cahn model. 


2.5 The jump of solutions of the transmission problems 

In this section we prove some results on the jumps of the solutions (n, T) and {u, T) 
of the transmission problems ()2.16p - (|2.20p and ()2.21l) - ()2.25l) . which we need in the 
following sections. 

We define a scalar product a :d f3 on hy a :d P = a : {Dj3), for a,P ^ S^. For a 
unit vector n G let a linear subspace of be given by 

<S^ = |-(a;(8 )n-l-n(8)a;)|cuGM^}, (2.84) 

let Pn '■ ^ be the projector onto 5)^, which is orthogonal with respect to the scalar 

product a :£) /3 and let Qn = I — P^- 

Lemma 2.10 Let uj* G be a vector. This vector satisfies ® u) — e)^n = 0 if 

and only if e{uj* 0 n) = Pn£ holds. 

This lemma is proved in [71 Lemma 2.2]. 

Lemma 2.11 Let {u,T) be a solution of the transmission problem (|2.16l) - (|2.20l) . As¬ 
sume that u is continuous in Q and that the limits exist and define continuous 

extensions ofS/xU from the set 7 ' to y'UT and from the set 7 to 7 Ur, respectively. Then 
we have 


[e(Vj,u)] = e{u* 0 n) 
[sPVj,!!)] — Pn^i 
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[f] = D{eiu* 0n)-e) 
[f] = -DQnS, 


(2.85) 

( 2 . 86 ) 
















Proof: Equation ()2.85p is proved in [ 8 l Lemma 2.2], (12.861) is proved in [7]. For com¬ 
pleteness we give the short proofs here. 

Since by assumption u is continuous across E and since VajU has continuous extensions 
from both sides of F onto F, the surface gradients and (Vru)(“) on both sides 

of F coincide, hence [Vru] = 0. Using the decomposition (12. 8 h and the dehnition (j2.85p 
of u* we therefore obtain 


[V^u] = + Vrh] = [{dnu) < 8 ) n] -|- [Vru] = [dnu] ® n = u* ® n. (2.87) 

Thus, by (|2.17D . 

Z)(e(u*0n)-e) =F»([e(V,u)]-e[S]) = [Z)(e(V,u) - e5)] = [f]. 

This proves ([2.850 . From (12.191) and (|2.85p we infer that 

0 = [T]n = (^D[e{u* <Sin) — e)^n, 

so that [e(Va;M)] = s{u* i 8 > n) = PnE, by Lemma ETOI Therefore we hnd 

[f] = D{[£{V^u)] -e) =D{Pne-e) = -DQ^e, 

which proves () 2 . 86 l) . ■ 

Lemma 2.12 Let {u,T) and {u,T) be solutions of the transmission problems ()2.16l) - 
(j2.20p and ()2.2ip - (I2.24p . respectively. Assume that u and u are continuous in Q and 
that the limits and {Vxu)^^'^ exist and define continuous extensions ofVxU and 

of VxU from the set 7' to 7' U F and from the set 7 to 7 U F, respectively. Then we have 


£ : T 

[dull] = [ m" 


[VxU] = [ 


£-.T 


r{s)^ 


u (Sn. 


-nsy 

Proof: The decomposition ()2.8I) yields 

[Vxu] = [(dnu) 0 n + Vru] = [{dnu) 0 n] + [Vrh] = [dnu] 0 n. 
From this equation and from (12.220 . (12.240 we infer 


( 2 . 88 ) 

(2.89) 


£:T 


£:T 


0=[f]n= i^D{[£iVxu)] - [^^]e)jn= (^T>(e([5„u] 0 n) - 

Thus, Lemma [2.101 the linearity of the projector P„ and (12.850 . (|2.860 imply 

e :T n * . 


n. 


£{[dnu\ 0 n) = [ f ' T jPne = [ f ' T ]e(u* 0 n) = e 


whence 


>"(5)- 


e : T 


"nsy 


risy 


\u 0 n 


-nsy 


dnu] - M*) 0 n -h n 0 ( [dnu] - 


e : T 


nsy 


= 0. 


We multiply this equation from the right with n and obtain 


([dnu] — + n(\dnu\ — ■ n = 0, 

V ^<(5) ' V"(<5) ' 

which means that [5„u] — ^ multiple of n. Scalar multiplication of the last 

equation with n yields ([5 „m] — [;^^]^*) ■ n = 0, from which we now conclude that the 
hrst equation in (|2.88p holds. The second equation is obtained from ([2.890 . ■ 
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3 The asymptotic solution 

This section forms the first part of the proof Theorem 12.31 We state in this section the 
form of the asymptotic solution The proof continues in Section 01 where 

we study properties of the functions So,...,S 2 appearing in the asymptotic solution. 
In Section [5] we use these properties to show that is an asymptotic 

solution by verifying that the estimates (|2.51l) " (|2.54p hold. This concludes the proof of 
Theorem 12.31 


3.1 Notations 

Before we can start with the construction of the asymptotic solution we must introduce 
more definitions and notations. In particular, we must introduce parallel manifolds to 
the manifold T and we must extend the definition of the surface gradients for functions 
defined on T, which are given in Section [2Tl to functions defined on the parallel manifolds. 
These definitions and notations are needed throughout the remaining sections. 

Let (5 > 0 be the number from (12.3p . For ^ satisfying —5 < ^ < 5 

r? = {{t,V + ^{1,1])^) I {t,r]) £ T} 

is a three dimensional parallel manifold of T embedded in Us, and 


r^(t) = {x £ ^ \ {t,x) £ T^} 

is a two-dimensional parallel surface of T{t) embedded in Us{t). Let ri, T 2 £ be two 
orthogonal unit vectors tangent to T^(t) at x £ r^(t). For functions w : r^(t) ^ M, 
W : r^(t) and W : r^(t) we dehne the surface gradient and the surface 

divergence on r^(t) by 


Vr^w = 

{driW)Tl -h {dr2w)T2, 

(3.1) 

Vr^IT = 

{dryW) (g) Tl -(- {dr2W) (g) T2, 

(3.2) 

divp^ W = 

2 

Tl .dr,W + T2-dr2W = Y,n- (Vr,IF)Ti, 

2 — 1 

(3.3) 

divr^ W = 

{dr,W)Tl + {dr2W)T2. 

(3.4) 


Clearly, with Vr dehned in (12.5p and (12.6h we have Vro = Vp. For brevity we write 

divr = divrQ. If w, W, W are defined on F^, we dehne : F^ Vr^IF : F^ i->- 

1^3x3, (iivrjVF : Fg 1 -^ M, divr^IT : F^ i-)> by applying the operators Vr^ and divp^ to 

the restrictions tci , W\ , W\ for every t. With these dehnitions we have the 
\r^(ty Ir^w’ Iryt) 

splittings 


Va;W{t,x) 

V.IT(t,x) 

diYxW{t,x) 


d^w{t, r], 0 n{t, T]) + Vr^w{t, r], ^), 
{dyW{t, r], 0) n{t, 7]) + divr^ W{t,ri,(), 


where, as usual, W{t,r],^) = W(t,r] + n{t, . 


(3.5) 

(3.6) 

(3.7) 
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The operators Vr and divp can be applied to functions defined on subsets of T. In 
contrast, the operator S/rj introduced next can be applied to functions defined on T x J, 
where J C R is an interval. For rciFxJ—^R, IF iFxJ—>■ R^, IF : F x J —> R^^^ 
consider the functions r] i-> j] i-)- IF 4 ^^(? 7 ) = W(t,r],^), r] i-)> = 

W{t,r],^), which are defined on F(t). To these functions the operators Vr and divr can 
be applied. We set 



= Vrwt,^{r]) € R^, 

(3.8) 

Vr,W{t,7J,C) 

= VrIFi,^(r/) € R='^^ 

(3.9) 

div^IF(t,77,0 

= divrIFt,.t(?7) € R, 

(3.10) 

div^IF(t,?7,0 

= divrII4,.c(?7) € R^. 

(3.11) 


If IF is defined on then {t,r],^) —)• = W{t,rj + n{t,r])^) is defined on 

F X (—5,(5). Consequently, the gradient V^VF is defined. The connection between V^VF 
and Vr,IF = VrpIFi is given by the chain rule, which yields 


V,IF(t, V, 0 = (Vre VF(t, rj + n(t, rj)^) (l + ^ V^n(t, t?)) . 


(3.12) 


In particular, we have V^IF(t, ? 7 ,0) = VrW{t,r]). Similar formulas and relations hold 
for div^IF, div^IF. If IF : R^ is constant on all the lines normal to F(t), for 

all t, we have W{t,ri,^) = W{t,ri). For such functions we sometimes interchangeably use 
the notations V,,IF and VpIF. Similarly, we interchangeably use the notations VrjW and 
Vpic, div^IF and divpIF, div^IF and divpIF if w and IF are independent of 
Note that by (13.121) we have for x € F^(t) that 

Vr,IF(t,x) = (V,IF(t,r/,0)A(i,r?,0, (3-13) 

where € R^^^ is the inverse of the linear mapping (/ + : R^ —>■ R^. 

From the mean value theorem we obtain the expansion 


= I + (3.14) 

where the remainder term RA{t, V: 0 ^ is bounded when (t, tj, ^) varies in Fx(-5,5). 
Insertion into (I3.13j) yields 

Vr,IF(t, x) = V^Wit, V, 0 {I + iRA{t, r], 0) . (3.15) 

For re : —>• R we consider Vr^u; and VrjW to be column vectors. For such w the 

equation corresponding to (I3.15P is 

Vr^w{t, x) = A^{t, r/, ^)S7rjw{t, rj,^) = {l + ^ R^{t, r?, ^)) V^rc(t, ij, (). (3.16) 

Furthermore, (13.3|) . (13.151) and (j3.inp together yield for IF : Zi^ —> R^ that 

2 

divr,IF = Y.Ti- ((V^IF)(/ + ^RA)Ti) = div^IF + edivr.jIF, (3.17) 

i=l 

with the remainder term 

2 2 

divr,glF(Z, r,,^) = J2n- ((V^IF)i?Av) = • {{VrWt,^)RATi) , (3.18) 

^=1 i=l 
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(3.19) 


and this equation implies for 1^ : —)■ that 

divr^H^(t, ry, = div^W + ^divr^^ll^, 

where divr,^!^ = j=ii^rjWt^^) Ti{Tj ■ RaTi)- The terms divr,.clT and divr,^lT are 
bounded when (t, ry,^) varies in T x (—<5, <5). 

For functions w with values in R we define the second gradients Vp^rc, by 

applying the operators Vr^, S/rj to the vector functions Vr^rc, S/rjW. For W with values 
in R3 we define second gradients Vp^bF, V^IF by applying these operators to the rows 
of Vr^bF, VrjbF. We remark that 


= divr^ Vr^ w 


is the surface Laplacian. 

Definition 3.1 Let I C M be an interval. For k,m G Nq and p = 1,3 we define the 
space 

C*=(/,C™(F,RP)) 

= {(t,?y,0 ^ : F X / —)> R^ I d^dlVjjW G (^(F x /), £ < k, i + j < m}. 

3.2 Construction of the asymptotic solution 

We start with the construction of the asymptotic solution We assume 

that the hypotheses of Theorem 12.31 are satisfied. In particular, we assume that there 
is a sufficiently smooth solution F = F^^^ of the evolution problem (12.141) . (I2.15j) . with 
so(T, Kr, si(u, T, T, 5o, 5i, defined in ()2.40l) - (I2.43F By this assumption, 

the function (u, T, u, T, Sq, Si) is known as a solution of the transmission-boundary value 
problem (I2.16p - (j2.32p . We use the notation 

l^(^) = ll*,’ l"(r) = 1 - l+(r), r^ = rl^(r). (3.20) 

10, r < 0 

Let fi G C^((—2, 2)) be a function satisfying 0 < (j){r) < 1 for all r G R and fiir) = 1 for 
|r| < 1. With the constant a from (j2.38p we define a function : Q ^ [0,1] by 

= 0, otherwise. 

By ()2.34p . cj)^x is equal to 1 in transits smoothly from 1 to 0 in vanishes 

in Qout ^. With this function the asymptotic solution is defined by 

U^^\t,x) = u[^\t,x) 4>^\{t,x) +U^ 2 ^\t,x) (1 - (j)i,\{t,x)), 

S'^^)(t,x) = s['"\t,x)(j)f,x{t,x) + S^'"\t,x) (l - (f>fj,x{t,x)), 

T^^\t,x) = D(^£{VxU^^\t,x)) —£S^^\t,x)'^, 
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(3.22) 

(3.23) 

(3.24) 



















where u'i \ are components of the inner expansion defined in lAs, 

and are components of the outer expansion [u^^\T 2 ^\ defined in Q \ F. 

The function is equal to the inner expansion in the 

region and equal to the outer expansion T 2 ^\ in the region Qo^'^- In 

the region both expansions are matched. 


The outer expansion The outer expansion is defined as follows. With the solutions 
(n, T) of the transmission problem (I2.16P - (|2.20l) and {u, T) of the transmission problem 
(|2.2ip - (|2.25p we set for (f, x) G Q \ F 


Up\t,x) 

= u{t,x) + p^‘^u{t,x) + ^u{t,x). 

(3.25) 

Sp\t,x) 

= S{t, x) + p-f‘^Si{t, x) + fiS 2 {t, x) + p‘f‘^S^{t, x). 

(3.26) 

Tp\t,x) 

= dU[V xVp'\t,x)) -eSp\t,xpj. 

(3.27) 


The functions u, Si,..., and another unknown function T solve the system of algebraic 
and partial differential equations 


-diVa;T = 0, 


f 

-f-.8 +^''{S)Si 
-f-.e + i;"{S)S2 + \^"'{S)Sl 

-f : e + ^"(5)S3 + ^'''{S)Si~S2 + 

b 


L>(e(V,n)-e52), 

0 , 

0, 


in the set Q \ F. 

u{t, x) 

r]) 


AV2 _ 

+- dtSi-\A,Si = 0, 

c 

Moreover, u satishes the boundary conditions 
= 0, {t,x) € [ti,t 2 ] X on, 

= \^/‘^u*{t,r]) f 5i(t,7?,C) - {t,r])£T, 

J—OO 'ip (1) 

= v) r Slit, v,C)- "'' dC 

Jo V ( 1 ) 

+ Aa*(Fr?) r ( f So{^)dd-e)dC, (t,r?)GF. 

J —OO ^ J —OO ^ 


(3.28) 

(3.29) 

(3.30) 

(3.31) 


(3.32) 

(3.33) 

(3.34) 


(3.35) 


Since by assumption F, T, T, Si are known from the evolution problem, this system can 
be solved recursively. To see this, note that (I3.30D yields 


= 


T:e 


ns) 

We insert this equation into (13.3111 and solve this equation for S 2 to obtain 

^ _ f:e p"{S) / r:£ 

^ P'{S) 2p'{S)^p'{S) 


(3.36) 


(3.37) 
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Using this function in (j.S.29p . we can determine u and T from the boundary value problem 
(I3.28p . (I3.29p . (I3.33P - (I3.35p . Finally, we can solve (I3.32p for S^. 


The inner expansion The inner expansion is essentially obtained by 

smoothing the jumps of the functions u, u and S from the evolution problem for the 
surface F. Before we can define the inner expansion we must therefore study in the next 
two lemmas the jumps of u and u across F. 

Let u* = [dnu] and a* = [d'^u] be the jumps of derivatives of u across F. These 
functions are introduced in (I2.35p . (I2.36h . For {t,x) = (t,x(t,r],^)) G Us we decompose 
u and it in the form 


u(t, x) 
u{t, x) 


+ v{t,x), 






ril) 


V^"(o) 


(3.38) 

(3.39) 


where ^ are defined in (j3.20p and where the remainder terms v and v are defined 
by (I3.38p . (I3.39p . The decomposition (j3.38p is motivated by the fact that 


Kv]=0, i = 0,1,2, 


(3.40) 


which follows immediately from (I3.38p . (|2.18p and (j2.35p . (12.361) . The first two terms on 
the right hand side of (j3.38p thus serve to separate off the jumps of the first and second 
derivatives of u at F. Similarly, the normal derivatives of first order of v do not jump 
across F. More precisely, we have the following result. 


Lemma 3.2 Let {u,T) and [u,T) be solutions of the transmission problems (I2.16P - 
(j2.20p and (j2.2ip - ()2.25p . respeetively. 

(i) V defined in (|3.39n satisfies 


[div] =0, f = 0, 1. (3.41) 

( a ) Assume that T is a -manifold. Suppose that u G 6*^(7 U 7', M^), it G 6*^(7 U 7', M^) 
and that u has -extensions, ii has -extensions from 7 to 7 U F and from 7' to 7' U F . 
With the function spaces introduced in Deftnition \3.1\ we then have 

V G C2((-<5,5),C'2(F))nC3((-5,0],C'i(F))nC3([0,5),C'Hr)), (3.42) 

i! G c^((-(5,(i),c2(F)) nc=^(F X (-^,0]) nc=^(r X [0,5)). (3.43) 


Proof: To prove (I3.4ip note that by definition of [dnw] in Section 12.11 and by definiton 
of in (I3.20|) we have 




F ; f (+) 

ra) 


(ff-) 


fi'fo) 


V V’"(i) V’"(o) ^ 


e : f - 


u 


From this equation, from the first equation in (|2.88l) and from (|3.39p we obtain (|3.4ip for 
f = 1. For i = 0 equation (13.411) is an immediate consequence of (I3.39P and ()2.23p . This 
proves (i). 
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(ii) Since F is a C^-manifold, the coordinate mapping !->■ {t,x{t,rj,^)) = 

(t, 7y+^ n(t, r/)) and the inverse mapping (t, x) i-;> {t,r]{t,x),^{t,x)) are It follows from 
this differentiability property of the coordinate mapping and from our differentiability 
assumptions for u that i-)- u{t,rj,^) is in F x (—(5,0] and in F x [0,(5), and 

that (t,r/) (->■ u*{t,'q) = n{t,r]) ■ [Vxu]{t,r]) belongs to C^(F) and {t,r]) (->■ a*{t,r]) = 
n{t,r]) ■ [dnVxu]it, rj) belongs to C'^(F). Since by (I3.38P we have 

v{t,V,0 = u{t,V,0 - 

these properties imply that 
1 

he Pi (c2+”^((-(5,o],c2-”^(f)) nC2+™([0,(5),C2-™(F))). 

m=0 

From this relation and from (|3.40p we conclude that (I3.42p holds. Relation (I3.43P is 
obtained in the same way, using (I3.4ip instead of (I3.40I1 . ■ 

For brevity in notation we define 

=s : De{Vxv{t,x)), a'{^) = d^a{t,r],^), (3.44) 

^(0 = <5'(t, 7],^) = e : De{Vxv{t, x)). (3.45) 

Later we need the following result, which shows how a{t, rj, 0) can be computed from the 
limit values of T and T at F. 


Lemma 3.3 The function a defined in (|3.45l) satisfies 


s: 

e : f 

(t(O) 


d(0)+e:[f] 

d(0)+e:[f] 


e ; f (+) 
^"(1) ’ 
e ; f(-) 
<(0) ’ 



(3.46) 

(3.47) 

(3.48) 


Proof: 

u*{t,r]) 


We apply the decomposition (12.81) of the gradient to the function W{t,r],^) = 


b"(i) 




e:f(- 


ip"{0) 


d- 


. This yields 


iVxW)^+\t,r]) = ((5nVF)(8)n + VrIF)^+^ =u*0n ^ I . 

V’"(l) 

(|3.39p thus implies 

* e:r(+) ^ 

(Vxti)^"^^ = u ®n — -h VxV. 

Insertion of these equations into (I2.22p yields 

_ , , £ : ^(+) 

= D{e{u* (8> n) — e) -h De{Vxv)- 

We take the scalar product with e on both sides of this equation and note (12.851) and the 
definition of d in (I3.45|) to obtain (j3.46p . Equation (I3.47|) is obtained in the same way. 
To prove (I3.48p . we add ^3.461) and (I3.47P and solve the resulting equation for (t( 0). This 
proves the lemma. ■ 
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Definition of the inner expansion We can now construct the inner expansion 
With the remainder terms v, v introduced in (13.3811 . (13.391) we set 
in the neighborhood Us of T 


+ v{t,x) + (3.49) 

S{^\t,x) = +^52(t,r/,^-^), (3.50) 

T^^\t,x) = D(^e{Vxu[^\t,x))-eS[^\t,x)J, (3.51) 

where by assumption Sq, Si are known from the evolution problem for F, and where the 
functions uq, ... ,U 2 are defined by 


«o(T??,C) 

= u*{t,r]) f So{'d)d'd, 

J —00 

(3.52) 


= u*{t,r]) [ Si{t,r],'d)d'd, 

Jo 

(3.53) 

U2{t,r],C,) 

= a*{t,ri) [ f So{'di)d'did'd, 

J —00 J —00 

(3.54) 


The function S 2 = S 2 {t,rjX) together with another unknown function si = si{t,7]) solve 
a boundary value problem. To state this boundary value problem let K{t, rj, ^) denote 
twice the mean curvature of the surface T^(t) at rj € rg(t). With the notation introduced 
in Section [22] we thus have At(t,r/, 0) = K,r{t,r]). We write At'(O) = d^K{t,r],0). 

The boundary value problem for S 2 and si consists of the ordinary differential equa¬ 
tion 

^"(*S'o(C))5'2(t,r/,C) - S2{t,r],C) = F2{t,V,C), (3-55) 

with the right hand side given by 


F2it,v,0 = <7(0) + e : [f] S, --e : {f)S[ - 

Cl Z 

+ XF^ (d'(O)C + e : De(a* ® n + Vru*) [ 5o(??) dA 

^ J —00 ' 

+ (^-Aac'(0)c)5', (3.56) 


and of boundary conditions. To formulate these boundary conditions, we choose y? G 
C°°(M, [0, Ij) such that 


set 


v^+(C) = 


<f(C) = < 

IV lA 

(3.57) 

<f(C) 

ray 

"-“Ano)' 

(3.58) 
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and define 


= v-(0 (e: f‘-> - »'(0)C 


2 V ^//(o) / 




r(i) 

+ ^ + Vrn*)C^ 


(3.59) 


With this function the boundary conditions are 

S2(t,7^,0) =0, (3.60) 

liin (S 2 (t,rf,(} - P 2 (t,ri,C)) = 0. (3.61) 

C^ltoo 

The function si = si(t, r]) in ()3.56l) is independent of C- It is determined in Section IT^ bv 
the procedure sketched at the end of Section which we apply of course to the boundary 

value problem (j3.55|l , (j3.56jl , (j3.60|) , ()3.6ip instead of the problem (12.2711 , (12.2911 - (12.3211 . 
The function si, whose explicit expression is given in (|2.41ll - (12.4311 . forms the second 
term in the definition (12.1511 of the evolution operator . 


4 The functions So,. .. , S 2 from the inner expansion 

The functions ,..., Sa in the outer expansion can be determined explicitly from (|3.30ll - 
(|3.32p , whereas the functions So,..., S 2 in the inner expansion are determined as solutions 
of three coupled boundary value problems to linear and nonlinear ordinary differential 
equations. It is not obvious that these solutions exist and what properties they have. We 
study these solutions in this section. 

4.1 The function S^o 

The first boundary value problem determining Sq is given by (I2.26p , (j2.28p , 

Lemma 4.1 Assume that the double well potential V’ satisfies (|2.38p . Then Sq is a 
solution of the boundary value problem ()2.26p . (j2.28p . if and only if Sq satisfies the initial 
value problem 

SUO = ^2^(So(()), CeR, So(0) = ^. (4.1) 

Proof: Let Sq be a solution of (I2.26p . (|2.28n . We multiply (|2.26p by Sq and obtain 

|(A«-I«)^)=o. 

or 

^(^ 0 ) - ^(5^)' = (4.2) 

By (j2.28|l we have lim^_).oo <S'o(C) = 1- From (j4.2ll and from (I2.38|l we thus obtain that 
lim^^co (5g(C))^ = —2(71. Using again (j2.28p . we infer from this limit relation that 
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lim^^oo >S'o(C) = Oj hence Ci = 0. We solve (I4.2jl for S'q and use that because of the 
boundary conditions (j2.28l) the function So must be increasing, hence S'q must be non¬ 
negative. This shows that a solution of (I2.26|] . (j2.28l] must satisfy the initial value problem 

To prove the converse we differentiate the differential equation in (jd.lji and obtain 
(I2.26P . We leave it to the reader to verify that the solution of (|4.1I) satisfies the boundary 
conditions (I2.28p . ■ 

Theorem 4.2 Assume that ip € 1],M) has the properties (I2.38p . Then there 

is a unique solution Sq € (0,1)) of the initial value problem (|4.ip . This solu¬ 

tion is strictly increasing and satisfies (j2.26p and (I2.28p . Moreover, there are constants 
Ki, ..., Ko > 0 such that for a > 0 defined in (12.3811 


0<5o(C)<i^ie-“l^l, 

— 00 < c < 0, 

(4.3) 

1 _ K 2 e-< < 5o(C) < 1, 

0 < C < 00 , 

(4.4) 

|5'5o(C)| 

— 00 < ( < 00 , i = 1,..., 4 . 

(4.5) 


This theorem follows immediately from the standard theory of ordinary differential equa¬ 
tions, and we omit the proof. 

Lemma 4.3 If ip satisfies the symmetry condition (j2.39p . then the solution Sq of (14. ip 
satisfies for all (p G M 


5o(-C) = 1-5o(C), S'q{C) = S'q{-C), (4.6) 

r-ICI 

/ So{d)d^= / So{^)dd + C^, (4.7) 

J —oo J —oo 

So{d)d'd-C'^ (4.8) 

J — OO ^ 

|^"(5o(C)) -^"(^(0)1 < (4.9) 


Proof: If the symmetry condition (j2.39l) holds and if Sq is a solution of the initial value 
problem (14.11) . then also C (l — 5'o(—C)) is a solution. To see this, note that (12.3911 
and (14.111 imply 

^2^(1 - Soi-O) = + (^ - ^o(-C))) = Soi-O)) 

= ^2^(5o(-C)) = (9c5o)(-C) = 9c(l - 5o(-C)), 

whence 1 — 5o(—C) satisfies the differential equation in (14. ip . Since we obviously have 
1 — 5 * 0 ( 0 ) = we see that 1 — 5o(—C) is a solution of (14. ip . Since the solution of this 
initial value problem is unique, we infer that 5o(C) = 1 — 5o(—C) holds, which implies 
s'oiC) = s'q{-0. 
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To prove (j4.7l) . note that (j4.6p implies for C > 0 


So{^) dd 


' —OO 

r-< 


So{d)di} + j 

SQ{d)d^+ r 
Jo 


So{d) d^ + 


So{'&)dd 

So{^) + So{-d)d^ 

So{^) + (1 - Soid)) dd 



So{d) dd + C- 


From this equation we immediately obtain (14.7p . The inequality (j4.8p follows from (I4.7p 
and from (14.3p . which yield 


So{d)dd-C^ 


-ICI 


So{d)dd 


»-ICI 


< 


Kie-^^^^dd = — 
a 


For the proof of (14.9p note that S{() = 1 for ^ > 0. Consequently, the mean value 
theorem and (|4.4I) together imply for C > 0 that 

|^"(5o(C)) -^"(5(0)1 = |^"(5o(C)) -^"(1)1 < |^"'(r*)(5o(C) - 1)1 < CiF 2 e-“l^l, 

with a suitable number r* between 5*0 (C) and 1. For (" < 0 an analogous estimate is 
obtained using (|4.3p and noting that 5(C) = 0 if C < 0. ■ 


4.2 The functions Si and S 2 


The solutions Si and S 2 of the second and third boundary value problems are studied in 
this section. The second problem determining Si is given by the equations (I2.27p . (12.2911 
- (I2.32p . the third problem, which determines S 2 , consists of the equations (13.550 . (|3.56p . 
(j3.60p . ()3.61ll . The properties of Si and S 2 , which we need in Section [5l are summarized 
in the next two theorems. 

To state the first theorem we need the function pi : F x M —R, which is defined by 
Pi{t,V,C) = £■ + e : (4.10) 


with p± introduced in ()3.58p . 


Theorem 4.4 Assume that belongs to (^^([0,1], R) and satisfies the assumptions (12.381) 
and the symmetry condition (I2.39p . Suppose that the function sq = so{t,r]) in (I2.32|) is 
given by ()2.40p . Let Sq be the solution of the boundary value problem (|2.26p . ()2.28p . which 
exists by Theorem \4.S\ 

Then for every {t,r]) G F there is a unique solution C Si{t,r],C) : R —>■ R of the 
boundary value problem (I2.27p . (12.290 ~ (I2.32p . The function Si belongs to the space 
C'^(R, (^^(r, R)). Moreover, there are constants Ki,..., K 3 such that for the constant a 
defined in (j2.38p and for all {t, ri,C) G F x R the estimates 


\D 




, < 

Ki, 

|a 

<2, 

(4.11) 

1 < 

K 2 e- 

-a|CI 

0 < j, |a| < 2, 

(4.12) 

< 

K^e- 

-a|CI 

l«l <2, j = l,2. 

(4.13) 


hold. 
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We do not give the proof of this theorem, since it is almost the same as the proof of 
Theorem 1.2 in [7]. Morover, it is obtained from the proof of the following theorem by 
simplification. The main difference between the two proofs is that the right hand side 
Fi of the differential equation (12.271) for Si is bounded, whereas the right hand side F 2 
of the differential equation (j3.55jl for S 2 grows linearly for ( —± 00 . 

Theorem 4.5 Assume that ip satisfies the assumptions given in the Theorem \4.4\ Sq 
be the solution of the boundary value problem (j2.26j) . (12.281) . and let Si be the solution of 
the boundary value problem ()2.27p . (12.291) - (12.3211 . Suppose that the function si = si{t,r]) 
in (j3.56p satisfies (j2.4ip with sio, sn given in (I2.42p . (I2.43p . 

(i) Then for every {t^rj) E T there is a unique solution ( 1 -^ S 2 {t,r]X) : M —M 0 / 
the boundary value problem (I3.55h . (j3.56p . (j3.6np . (I3.61h . The function S 2 belongs to 
C^(M, R)), and there are constants K ^,..., Kq such that for the constant a defined 

in (j2.38p and for all {t, ry, C) E T x R the estimates 

\diDl^^^-^S2{t,r],C)\ < i^4(l + |CI)'"', |a|<2, j = 0,l, (4.14) 

l^^h, 77 )('S' 2 (T?/,C) - P 2 (T??,C))| < -^ 5(1 + 0<j,\a\<2, (4.15) 

\dcD'^t,r>)S2{t,7],C)\ < i^ 6 (l + |CI)e-“'^', l«l<2, (4.16) 

hold, where p 2 is defined in (I3.59h . 

(a) S 2 is the only solution of the differential equation (j3.55p with F 2 given by (j3.56p . 
which satisfies (I3.60h and for which constants C ,6 > 0 exist such that 

|S 2 (t,r?,C)l<Ce(“-^)l^l, CGM, (4.17) 

holds. 

4.3 Proof of Theorem 4.5 

In this section we give the proof of Theorem 14.51 which is divided into five parts: 

(I) Reduction of the boundary value problem for S 2 to a problem in L^. With 


P 2 defined in (j3.59p we make the ansatz 

S 2 {t, r], C) = h, C) + P 2 {t, V, C)- (4.18) 

Insertion of this ansatz into the equations (I3.55P and (I3.60p . (|3.61l) shows that S 2 is a 
solution of the problem given by these equations if and only if w solves the equations 

i’”{So{C))'w{t,rjX) - d‘^w{t,ri,C) = F 2 {t,p,C) + F 3 {t,ri,C), (4.19) 

w{t,p,0)=0, (4.20) 

lim w{t,r],C) = 0, (4.21) 

C—>■±00 

where F 2 is given by (j3.56p and where 

Fs = -{^"{So) - df)p 2 . (4.22) 


To get (I4.20|) we used that <y 2 +( 0 ) = v^-(O) = 0, which by (I3.59p implies p 2 {t,ri, 0 ) = 0. 
To show that the solution S 2 of the problem ()3.55l) and ()3.60p . (I3.6ip exists, it therefore 
suffices to prove that the reduced problem (14.191) - (I4.22p has a solution. 
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(II) Spectral theory For this proof note that is a linear self-adjoint 

differential operator in L^(]R). From the spectral theory of such operators we know that 
the continuous spectrum of V’"(5o) — is contained in the interval [oq, oo), where 

ao = min| lim ^"(5o(C)), hm V;"(S’o(C))|, 

and that the part of the spectrum in (— 00 , oq) is a pure point spectrum. (14.9h yields 
lim '0"(S'o(C)) = lim i^"(S'o(C)) = hence the assumption ()2.38p implies 

C—>■ —00 <^—>-co 

ao = > 0. Therefore 0 does not belong to the continuous spectrum. From the spectral 

theory we also know that for every w G C, which is not in the continuous spectrum, the 
differential equation (i/;"(5o) — — uw = / has a solution w G if and only 

if / G L^(M) is orthogonal to the kernel of the operator 'ip"{So) — d‘^ — u. This implies 
in particular, that for every {t, 77 ) G F the differential equation (14.191) has a solution 
w{t,r],-) G T^(M), if the right hand side C /(C) = -^ 2 (^ 5 ^jC) + belongs to 

L^(M) and is orthogonal to the kernel of the operator To show that the 

problem (j4.19p - (j4.2ip has a solution, we therefore verify in the next two parts of the 
proof that F 2 -|- T 3 satisfies these two conditions. 


(Ill) The asymptotic behavior of F 2 + F 3 at infinity. We first show that the right 
hand side F 2 + F 3 of (I4.19p decays exponentially at ± 00 , which implies that F 2 -|- F 3 G 
L^(M). To simplify the notation we define 


^+(0 = ^+(C)^"(5o(C)), ^-(C) = v:>-(C)<(5o(C)), (4.23) 

with (/?_!_, (f- given in ()3.58p . For these functions we obtain from (|4.9p that 

|^_ - ^"(0)(^_| = <^_(C) |^"(5o) - (4.24) 

= ^+(0 \X'{So) - ^"(1)1 < CiF4e-“l^l, (4.25) 


for all C € M. Since tp "= 1 for C < —2 and V’"(1 )(/j+(C) = 1 for > 2, these 
estimates imply 

1 1 - (01 < CiF4e"“l^l, -00 < C < 0. (4.26) 

|1-^+(01 < CiF4e"“l‘^l, 0<C<oo, (4.27) 

|l-^_(C)-^+(C)l CgM. (4.28) 


To get the last estimate we combined the first two estimates and noted that ^_(C) = 0 
for C > “1 and ^_|_(C) = 0 for ^ < 1 . 

Note that by (I3.56p . (j3.59l) and (14.221) the function F 2 + F 3 can be decomposed in the 
form 


F2 + F3 


F2-^"{So)p2 + d^^P2 = '£^F 




(4.29) 
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where 


2 ^ 2 V w/hn'i / ^ 2 V ^w(i) ; ^+’ '' ^ 


Ii = <T(0)+e: [f]5i-e:f(-)^_-e:f(+)^+, 

h 
h 
h 


(4.30) 


V’"(0) 

<t'(0)C (1 -^_ -^+), 


— 


e : De{a* (g) n + VxU*) j 5'o(i?) dr? — i 

— — e : {T)S[ + — Ak'(0)C^ Sq + 9^/32 • 


(4.32) 

(4.33) 

(4.34) 


We show that everyone of these terms decays to zero for () —)• ±oo. To verify this for the 
first term we insert (j3.46p and (j3.47p into (I4.3n|] . which results in 


I, = a(0) + e : [T] 5i - <t(0 )(^+ + ^-) - e : [T] 


F : f (-) 

^"( 0 ) 


+ 


F ; f ( + ) 

^"( 1 ) 




We introduce the terms and into this equation. Noting the definition 

of Pi in (I4.10p . this leads to 

l^il < |'3-(0) {l-Tp^+Tp_) + e: [f] (S’! - pi)\ 


+ 


e:[T] 


e : f(-) 


{(f- --^"( 0 )(/?_) + 


e : f (+) 


(<7’+ - i’’' 


<C'e““l^l, (4.35) 


V>"( 0 ) V'"(l) 

for all C £ 1^) where we applied the estimates (I4.12p . (14.241) . (I4.25P and (|4.28l) . 

Next we estimate l 2 - By definition we have ^+(C) = 0 for C ^ 1- From (I4.3ip we 
thus have on the half axis — oo < C < 0 that 


h = - 




+ 


■I 


V’"(o) 




r{Q)-r{So)^2 r{^).2r. <'(o) 

=-^--Si - - -e-) - 

= hi + h2 + -^23- 


sl- 


e : 


V;"( 0 ) 




(4.36) 


To estimate /21 we apply the mean value theorem to i!)'" and use (j4.3p and (j4.1ip . to 
estimate h 2 we use (14.261) and (14.lip . The result is 

1/21 + / 22 I < CKie-^\^\ -00 < C < 0. (4.37) 

To estimate h'i note that by (j3.57p . p3.58p and (I4.10p we have for —00 < C, < —2 that 

e : fh)[t,p) 


Pi{t,pX) = ■T'' \t,p)) = 


V’"(0) 


With this equation we infer from (14.111) and ()4.12p with a, j = 0 that 


1^231 = 


V^"'(o) 


5i- 


e : fh) 
^"( 0 ) 


+ -cx)<C<0. (4.38) 

7//Tn'i J 


V^"( 0 ) 
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(j4.36p - (I4.38|) together imply that |/ 2 (C)I < Ce “1^1 for —oo < C < 0- On the half axis 
0 < C < oo we estimate I 2 analogously. This proves that 

I/ 2 I < —00 < C < 00 . (4.39) 

The estimate for is obtained by application of (j4.28j) to (j4.32j) . which immediately 
yields 

l^fsl < 0(1 + -00 < C < 00 . (4.40) 

To study the asymptotic behavior of note that (14.81) and (I4.27p together imply 


5o(T9)di?-C+^+(C) < r So{^)d^-C+ +c+|l-^+(C)| 

' j —OO 


< [-Ki + CCKi]e 


-a|CI 


Insertion of this inequality into (I4.33h results in 

|/4|<0(l + |C|)e-“l^l. 


(4.41) 


It remains to investigate Note first that the third term on the right hand side of 
(j4.34p satisfies 

dlp{t: V, C) = 0, for Id > 2 . (4.42) 

To show this it suffices to remark that the functions (p± are constant on the intervals 
(— 00 , —2) and (2, 00 ), from which we see by inspection of (j3.59p that on these intervals 
the function ( 1 — p 2 {t,r],Q is a sum of constant and linear terms, whence (|4.42l) follows. 
If we estimate the first term on the right hand side of (I4.34p by employing (I4.13j) with 
a = 0, j = ^ and the second term by using ()4.5I) . we obtain together with (I4.42p that 

I/ 5 I <C(l + |d)e-“l^l. (4.43) 

We combine ()4.29l) . ()4.35p . (14.3911 - ()4.41l) and (|4.43p to derive the estimate 

|F 2 (Tr?,C)+T 3 (Tr?,C)| <C(l + |C|)e-“l^l, for all C G M, (4.44) 

which shows in particular that the right hand side of (j4.19p belongs to L^(M). 


(IV) The orthogonality condition determining si. Next we must show that the 
right hand side of (|4.19l) is orthogonal to the kernel of — This kernel is different 

from {0}, since Sq belongs to the kernel. This is immediately seen by differentiation of 
(I2.26p . which yields 

i^"{So{0)S',{^)-d^S',{C)=0. (4.45) 

Since by (14.5p the function Sq is in the domain of definition of 'ip"{So) — d^, it belongs to 
the kernel of this operator. 

The theory of linear ordinary differential equations of second order implies now that 
the kernel is one-dimensional, hence every function from the kernel is a multiple of Sq. 
Therefore the right hand side F 2 +F 3 of (|4.19l) is orthogonal to the kernel if it is orthogonal 
to Sq. Note that the integrals -^ 2 <S'o FqSqcIC both exist, since F 2 and T 3 
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grow at most linearly for ^ —>■ ±oo, whereas by (14.51) the function Sq decays exponentially 
at ±oo. Therefore we obtain from (j4.22p and (I4.45P by partial integration that 


/ OO /*00 /* OO 

(F 2 + F3)5' dC = / F 2 S'odC- (- df)P 2 ) S'o dC 

-OO j —OO j —OO ^ ' 

/ OO poo poo 

F2S'odC- p2{i>"iSo)-dl)S',dC= F2S',dC. (4.46) 

-OO J —OO J —OO 

To study the last integral on the right hand side note that by (14.6h the function Sq is 
even, which implies that 

/ OO poo 

(T'(t,77,0)C5’o(C)dC = 0, and / ^'(t, t/, 0)C (5’o(C))^dC = 0. (4.47) 

-OO J —OO 

Moreover, since by Lemma l4.1l the function Sq satisfies (14.ip . we obtain by substitution 
of ^ = 5o(C) 

/ OO poo j p\ j 

S'Q{OS'Q{C)dC= J2^{SoiC))S'Q{OdC= V2^(7?)dr? = ci, (4.48) 

-OO J —OO J 0 

where the last equality sign holds by definition of ci in (|2.37p . Finally, by partial inte¬ 
gration, 

r So{d)ddS'Q{OdQ= \\xn ( So{d)ddSo{(:i)- 5o(C)"c) 

J —OO J —OO J —OO J —OO ^ 

= hm SQ{C){So{C,)-SQ{0)dC 

/ OO poo 

So{C){l-So{C))dC= SoiC)So{-C)dC. (4.49) 

-OO J —OO 

In the second last step we used that Sq is increasing. The equality sign thus follows from 
the theorem of Beppo Levi. The last equality sign is obtained from (14.6p . 

The equations (|3.56l) and (|4.47l) - (|4.49l) yield 

/ OO poo 1 poo 

F2{t,p,C)S'o{C)dC= / {a{0)+e:[f]Si)S'QdC--e:{f) / S[S'QdC 

-OO J —OO J —OO 

1 ^00 poo p(^ 

- ^ / i’'”iSo)SfS’Q dC + A^/2 £ . 0 n + Vr^*) / / So{^) d^ S'q{C) dC 

^ J —OO J —OO J —OO 

+ r (AVV(0)C + (--A«'(0)C)5(,)5(,(C)dC 

J—OO ^ C / 

/ OO 1 poo 1 poo 

S,S'q dC--e: (f) / S[S'q dC-- ^"'(5o)5?5(, dC 

-OO J —OO ^ J —OO 

*5o(C)*5o(-C)dC + ci-. 

-OO ^ 

= Sio - A^/^—Sii + —Si = — (-S 10 - A^/^sii + Si). 
c c c c 

To get the second last equality sign we inserted (|3.48p for (t( 0) and used ()2.42l) . ()2.43l) . 
The right hand side of this equation vanishes if and only if si satisfies (I2.4ip . From (I4.46p 
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we thus infer that F 2 + -F 3 is orthogonal to the kernel of -i/;"(S'o) — if and only if (|2.4ip 
holds. 

Consequently, from part (II) of the proof we conclude that the differential equation 
(I4.19P has a solution w in L^(M) if and only if si satisfies (|2.4ip with sio and sn given in 
(j2.42p . p.d.ip . In fact, there is exactly one such w, which also satisfies (j4.20p . To prove 
this assume that w & L‘^ (M) is a special solution of (I4.19P . Then we obtain every solution 
contained in L^(M) in the form w = w + (33^ with an arbitrary constant /3 G M. Since 
(j4.ip and (I2.38P yield 

5 '(0) = 72^(50(0)) = ^2^(^) > 0, 

we can choose /? such that 


w{t, r], 0) = w{t, 77 , 0) + /35o(0) = 0, 

which is (I4.20p . This equation determines f3 uniquely, hence w is the unique solution of 
(j4.19p and (j4.20p in L^(M). 

(V) Existence of the solution, estimates (4.18) — (4.20). We show next that this 
function w satisfies (I4.21I) . To this end we need the following 

Lemma 4.6 Let a_ > 0, a+ > 0 and set d = ^ min{a_, a+}. Let (7 : M — M 6e a smooth 
function and /et / : M — > M 6e a continuous function such that 


b(C)-a-| < 

for C < 0, 

(4.50) 

b(C)-a+| < 

for C > 0, 

(4.51) 

1/(01 < C(l + |C|)e-“l^l, 

for C G M. 

(4.52) 

Let w be a solution of 

g{0w{C) -d^w{C) = /(C), 

c G M. 

(4.53) 


(i) Then w belongs to the space C'^(M). If w ^ L^(R), then there is C > Q such that 

|i9^Ti(C)l < ^*(1 + /or all C £ M, /or j = 0,1, 2. (4.54) 

(a) If there are C,6 > 0 such that 

|h;(C)| < (4.55) 

holds for all C G then w € L^(M). 

This is a standard result from the theory of ordinary differential equations, and we omit 
the proof. 

To show that w satisfies (j4.2ip . we apply this lemma with a_ = '!/)"(0), a+ = 
^"(1), 5(0 = '0"(5'o(C)) and /(C) = F 2 (t,r],C) + F3{t,r],C)- Then we have o = 

min{Vj"(0), 7^"(1)} = a, by (12.381) . and from (14.91) and (14.441) we see that (I4.50p - 
(I4.52P hold for this choice of functions and constants. Moreover, with this choice of func¬ 
tions the differential equation ()4.53p is equal to (j4.19p . Since w G L^(E) is a solution of 
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(j4.19p . we see that all assumptions for part (i) of Lemma 14.61 are satisfied, hence (j4.54p 
holds for w, which means that 

|5'^u;(C)| <C(l + |C|)e-“lfl, CGK,j = 0,1,2, (4.56) 

and this in particular implies that w satisfies (|4.2ip . 

We have now found a unique solution w G L^(M) of (I4.19h - (I4.21I) . By part (I) of 
this proof this means that S 2 given by (14.181) is a solution of the boundary value problem 
(|3.55p . (|3.56p . (|3.6np and (I3.61h . Since by (j4.18p we have w = S 2 — P 2 , the inequality 
(I4.56P shows that S 2 satisfies ()4.15p for a = 0. 

To verify that S 2 satisfies (I4.15p for a / 0, it must first be shown that S 2 is two 
times continuously differentiable with respect to This follows if we can show that 

w = S 2 — P 2 two times continuoulsly differentiable with respect to (f,r/), since by our 
regularity assumptions the function p 2 has this differentiability property. To prove this 
differentiability of w, we write (I4.19p . (|4.20p as a perturbation problem for the linear 
equation 

Aw = f{t,r]) 

in L^(M), where A = {'ip"[So) — d^) is the linear differential operator on the left hand side 
of (j4.19l) and fit,r]) = F 2 {t,r], •) + F 3 {t,rj, ■) G L^(M) is the function on the right hand 
side of (I4.19p . which depends two times continuously differentiable on (t,r]) and satisfies 
the estimate (j4.44l) for every G T. Since 0 is an eigenvalue of A and since f{t,rj) 

is orthogonal to the kernel of A for every this linear equation has infinitely many 

solutions and the solution set is affine. The condition (|4.20l) defines a linear subspace, 
which is closed in the Sobolev space and which intersects the solution set in 

exactly one point w, which is the solution of (I4.19p . (I4.20p . 

To the problem set in this way we can apply the pertubation theory of linear operators. 
The theory yields that w is two times continuously differentiable with respect to (f, ??). 
We avoid the details but refer to standard texts on the pertubation theory of linear 
operators, for example [26] . 

With this knowledge we can derive the estimate ()4.15p for a 7 ^ 0 by applying the 
differential operator with 1 < |q;| < 2 to the differential equation (I4.19P and obtain 

riSo){Dl^)w) - = DI^^{F2 + F,). (4.57) 

This is a differential equation for the function w with right hand side satisfying the 
estimate 

|^ft,,)(i "2 + ^3)1 < C{1 + |C|)e-“l^l. (4.58) 

The proof of this estimate proceeds in the same way as the proof of the corresponding 
estimate for a = 0 , which we gave in part (III). Essentially one has to replace the terms 
appearing in F 2 and T 3 , which depend on by their derivatives. To avoid repetition 

of many technical details, we omit this proof. 

The differential equation (I4.57P has the same form as the differential equation (|4.19l) . 
From (I4.58P we see that the assumption (|4.52p holds, hence we can apply Lemma [4.6l (i) to 
this differential equation, from which we see that w belongs to the space (7^(1^, C'^(r,R)) 
and that the inequality (I4.54p holds with w replaced by whence we have 

\diDl^-^w{C)\ < C(1 + |C|)e-“lfl, C G M, 0 < i < 2, 1 < |a| < 2. 
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Since w = S 2 — P 2 , this is inequality (I4.15p with a 7 ^ 0. Therefore we proved that (j4.15p 
holds for all 0 < |q:| < 2 . 

The inequality (14.1611 is a consequence of (14.1511 and of (14.4211 . the inequality (I4.14p 
follows by combination of (14.1511 with the estimate 

<C(l + |C|)^■^ i = o,i, 

which is seen to hold by inspection of (I3.59p . 

We have now proved statement (i) of Theorem [431 ^.nd it remains to verify (ii). That 
is, we have to show that S 2 is the only solution of p.55p . p.56p . ()3.60p satisfying ()4.17ll . 
Indeed, from (I4.14h it follows that S 2 satisfies (j4.17p . Assume that S'! is a second solution 
satisfying ()4.17ll . Then w = S 2 — S 2 fulfills (|4.55ll and the equation 

— d'^w = 0 . 

Lemma 136] fiil thus yields w € L^(]R). Consequently, by (|4.18p we have = w + w + p 2 -, 
where w + w ^ L^(M) is a solution of (j4.19p . (I4.20p . At the end of part (IV) of this proof 
we showed that w is the only solution of (|4.19ll . (I4.20p in whence ib = 0, hence 

= 52. 

The proof of Theorem 14.51 is complete. ■ 

5 Proof of the estimates (12.511) — (12.541) in Theorem 12.31 

The proof of (12.5111 - (I2.54p is straightforward: We insert the function 5^^^) 

defined in Section [3| into the model equations (jl.ljl and (jl.3p and compute the residues. 
However, the necessary computations are long. Therefore we divide them into four parts: 

In Section [5.11 we compute for the functions and from the inner and outer 
expansions of the residues diva;T^^^^ + b and + b separately. Likewise, 

in Section 15.21 we insert the inner expansion s[^^) and the outer expansion 

into don and compute and estimate the residues separately. With 
these residues we can prove (|2.5ip - ()2.54p in the regions and but in the 

matching region we need auxiliary estimates, which are stated in Section [531 All 

the estimates are put together in Section 15.41 to complete the proof. 

5.1 Asymptotic expansion of divajT^'^) + b 

Lemma 5.1 Let {u,T) be the solution of the transmission problem (I2.16h - (I2.20p . With 
the splitting (I3.38P of u the stress tensor field T satisfies in the neighborhood of T 

f = [f]5 + DeiVxv) + De{a* (g) n + Vr^u*)t + DeiVr^a*)^{^+f. (5.1) 

With a' defined in (|3.44l) and with a remainder term R-.f € L°°{Us) we have for ft, 

Us 

e : T{t,P,C) = e : + a'{t,p,0)f 

+ £:De{a*{t,r])^n{t,p) + Vru*{t,p))f, + R^.f{t,p,(,)f,‘^, f, > 0, (5.2) 

£-T{t,P,C)='S-T^~Ht,r]) + a'{t,r],0)^ + R^.f{t,p,C)f, f < 0. (5.3) 
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The functions u* and a* introduced in (j2.35j> . (j2.36p and the normal vector n satisfy on 
the interface T 

(^D[e{u* ^ n)— e))'^n = 0, (5.4) 

[De{a* 0 n+ Vru*))n + divrDe{u* (S'n) = 0. (5.5) 

Proof: With the splitting (13.611 of the gradient operator we compute from (13.381) that 

Vxu{t, x) = d^u (g) n + Vr^n 

= {u*{t, ??)l+(0 + ® n{t, rj) + Vr^u*(t, r/)^+ 

+ Vr^a*{t, + 'V^vit, x). 

We insert this equation into (12.171) . note that S{t,x) = l'''(C)) employ that by (12.851) 

D{eiu* (g)n)-e)S = [f]5 

to obtain (j5.1|) . 

By (13.151) we have 

Vr,u*{t,v)t = Vru*{t,v)f^ + (V^u*(^,7?))i^A(^,^?,0(?+)^ (5-6) 

where we used our convention to identify Vru* and V^n*, since u* does not depend on 
Noting the definition of a in (13.441) . we obtain from (|5.ip and from (15.61) that 

e : = e : [T]{t,v)'^^iC) + 

+ e : De{a*{t, rf) ® n{t, r?) + Vru*(t, r/))^+ 

+ e : De{Vru*{t, r])RA{t, rj, ^ + ivr 5 a*(t, rj)) (^+)2. (5.7) 

(j3.38p and (j2.17p together imply for .^ < 0 that 

T{t,r],C) = De{yxv{t,r],^)), (5.8) 

whence 

e : T{t,r],^) =e : De{\/xv{t,ri,^)) = &{t,r],^), ^ < 0, 

and therefore 

s :f^~\t,r]) = a{t,r],0), s : f^+\t,r]) = s : [f]{t,r]) + a{t,r],0). (5.9) 

By Taylor’s formula we can express a in the form 

o-(t, V, C) = r/, 0) + d^a{t, ry, 0)^ + 9|d(t, ry, C)^^- 

We expand a in (15.7p with this formula and note the equations (15.9p to obtain (15.21) and 

dES]). 
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()5.4p is an immediate consequence of (12.851) and (12.191) . To prove (j5.5l) we apply (13.7p 
to calculate from (15.ip that 

0 = div^T + b 

= d^{[f]nS) + dWr^[f]S 

+ d^{De{a* <^n + VrjU*)^+)n + d^[^De(yr^a*){^'^)^)n 

+ divr,Del^{a* ® n + Vr,u*)t + ^(Vr,a*)(e+)2) 

+ div^DeCVxv) + b. (5.10) 

From this equation we obtain for ^ < 0 that 

diVxDeiVxv) + b = div^T + b = 0. (5-11) 

By assumption in Theorem 12.31 the function b is continuous at T. Moreover, by 
Lemma 13.21 and the differentiability properties of v required in Theorem 12.31 the function 
V is two times continuously differentiable at T. Therefore we infer from (15.lip that 


(diVajHepVa;?)) + b)^"'"^ = [diVxDe{Vxv) + b)'' "’ = 0, on T. 

With this equation and with [T]n = 0, by (I2.19h . we conclude from (I5.10h that 
0 = ^lim (div^T + b) = (^De{a* (g) n + Vru*))n + divr[T]. 

From this relation and from divrfT] = divrLle(u* 0 n), which is a consequence of p2.85p . 
we obtain (15.5p . ■ 

Next we study the stress field in the inner expansion. 

Lemma 5.2 Let u[^\ s[^\ he given in p3.49p - p3.5ip . let uq, ui, U 2 be defined in 
(I3.52P - (I3.54p . and let Ra be the remainder term from (13.1411 . We set ( = Then 

we have for {t,rj,f^) from the neighborhood Us o/F 

T[^\hg,0 = [f][So + gi^/^Si)+De[Vx{v + gi^/^v)) 

+ ® n + Vru*) [ Soiid) dd + (5.13) 


(-) 


(5.12) 


where 


i?Ti(A,Lr?,^,C) = ^(e(Vre(A^/'ui + Au2) + AC(V^uo)i?A) -£^ 2 ). (5.14) 

The argument o/[T], u*, a*, n is {t,r]), the argument of Si, S2, uq, ui, U2 is {t,r], (^>,^1/2 ); 
the argument ofVxV, ^xV, Ra is {t,rj,f), and the argument of Sq outside of the integral 


is 


Moreover, we have 


(/iA)i/2 • 

divj; t[^'’ + b = diVj,L>e(Va;i)) + b + ^ divr,5[T]5'o + /r^/^divr^ ([^’I'S'i) 

+ ti^^‘^divxDe{Vxv) + iiiX)^^‘^divr^De{a* ^ n+ Vru*) f So{d)dd 

J —00 

+ g,dYVxRTi- (5.15) 
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With o'(O), o''(0), d'(O) defined in (13.441) . (j3.45p we have 


= -e : = -e : + - &{0) - {0)C - //'/^(O) 

— : De{a* ^n + Vru*) ( So{'&)d'd — fj,RwiX,t,r},^X): (5.16) 

J —OO 

where 


Rw{X, t, ?7, C) = e : Rti (A, t, rj, C) + Xd^d{t, r], + X^/‘^d^d-{t, r], |)C, (5.17) 

with suitable ^ between 0 and ^ and with Rt^ defined in ()5.14p . 

Proof: With the splitting (j3.6|) of the gradient operator we obtain by definition of 
in (IMP . (IMP - (IMP that 

= ifJ-X)^^'^VxUo +tiX^^^VxUi+fiXVxU2 + Vxiv + 

= {u* ^n)So +^nSi +{iaX)^^‘^a* f So{'&)d'd 

J —OO 

+ (/uA)^/^Vrjtio + ^A^/^Vr^ui + /uAVr^n2 + Vx{v + (5.18) 

(|3.52p and ^3.151) together yield 


Vr,uo = V^no(/ + ^Ra) = Vpn* So{^)d^ + (^A)^/^^ (V^uo)i?A. (5.19) 

J —OO 

We insert (|5.18P . (I5.19P and (|3.50l) into (|3.51P . From the resulting equation we obtain 
P5.13p and ^5.141) if we also note that by (12.851) 

D[e{u* ® n)(5o + - e{So + 

= D{e{u* ®n) -e){So + = [f]{So + 

To prove P5.15P we employ the splitting (13.7p of the divergence operator and (I3.19P to 
compute from (I5.13P 

div^T,^'^) + b = a^([f]n5o) + (divr[f])5o + e(divr,g[T])5o 
+ ^x^/^{^^{[f\nSl)+divr,mSl)) 

+ diWxDeiyxv) + b + ^^/“^dAVxD^iyxv) + [De{a* ® n + Vru*))nS'o 
+ (/iA)^'^^divrjT)£(a* 0 n + Vru*) f So{'d)d'd + fvdivxRTi- (5.20) 

J —OO 

By (|2.85l) and (15.51) we have 

(Zle(a* 0n +Vru*))nS'o + divr[r]5'o = ^(lle(a*0n + Vru*))n + divrIle(u*0n)^S'o = 0. 
With this equation and with p2.19P we obtain p5.15P from p5.20p . 
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(|5.16p . (|5.17p follow immediately from ()1.8I) . which implies dsy\l{£{'Vxu[^^),s[^'^) = 
—e : , and from ()5.13p , ()5.14p , using the Taylor expansions 

(t(0 = <t(0) + d^cr{i)i = ,t( 0) + (/xA)^/29^d-(|)C, 

(t( 0 = <7(0) + 5go-(0)^ + dla[i)f = (t(0) + (0)C + ^A(t"(|)C^. 

This completes the proof of Lemma 15.21 ■ 

Corollary 5.3 Let ^ and defined in (|2.34l) and let he given by (I3.5ip . 

Then there is a constant C such that for all 0 < fi < fio and all 0 < X < Xq 

1 /^2 

Proof: We estimate the terms on the right hand side of p5.15p . Note first that if 
(L ri, 0 G U and C = then 

1^1 <-(/iA)^/^|ln/r|, |C|<-|ln^|. (5.22) 

a a 

This follows from ^2.341) . With these inequalities the hrst two terms on the right hand 
side of (I5.15P can be estimated as follows: From the differentiability properties of v and 
b, which in Theorem 12.31 are assumed to hold, it follows by Lemma 13.21 that 

div^L>e(V,,i;) + b G C{Us) n ((-<5, 0], C(r)) n ([0, 6), C(r)) . (5.23) 

Because of this differentiability property we can apply the mean value theorem to 
diVxDe{S/xv) + b, which together with (|5.12l) and (|5.22l) yields for all (t, G U 
Qm^ch with e > 0 that 


{dwxDeiVxv) + b) {t, rj, ^ = (divo; De{Vxv) + b) +5^ {d\NxDe{Vxv) + b) (t, r], C)f, 

= \d^{diVxDe{Vxv) + b)(t,r/,^*)|^ < < C'i^(/xA)^/^| ln/r|, (5.24) 

with a suitable number between 0 and Since by (15.111) the term diYxDeCVxv) + b 

vanishes for C < Oj the inequality (I5.24h holds for all (t, ?7,0 ^ Qinn^ U Q^atch- 
To estimate the last term in (15.151) note that (14. Sp . (14. lip and (|5.22l) yield 




0 < 


' —oo J —oo 

f 

lo 


0< ^ So{d)dd<C++ C2<C3-\lnn\, 

J — OO ^ 

So{di)dididd < -(C+)" + C4 < C5(-|ln/i|)^ 


(5.25) 


Si{t,r],'d)d'd < Ce\C\ < CG-\ln fi\. 


Using these inequalities, the definitions of uq, ui, U2 in p3.52p - (I3.54p and the inequality 
p4.14p we obtain from (I5.14p that 

Cs. 


\RtA < (C7 + ^)|ln^|2, 

a2 


|^diVa:i?Ti| < A (Cg + ^) | In/ip 
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(5.26) 

(5.27) 

















































(j5.26p is used later, (|5.27|) is the desired estimate for the last term in (j5.15p . 

To estimate the other terms in (j5.15l) we apply (|4.3I) . (|4.4I) . (|4.1ip . (15.221) and ()5.25p . 
Together with (I5.24h and (j5.27p we find for G U Q^atch 


(diva^r^^^^ + h){t,r],^) 


< 




+ 


<^12) + (AiA)^/^|ln/r| + ^/^)|ln/rp 


which implies (|5.21l) . ■ 

In the next lemma we study the outer expansion . 


Lemma 5.4 Let {u, T) he the solution of the boundary value problem, which consists of 
the elliptic system (I3.28p . ()3.29p with S 2 given by (13.371) . and of the boundary conditions 
P3.33p “ p3.35p . Let S 3 be the solution of p3.32l) . Then defined in p3.27p satisfies on 

Q\r 


= f + + plT - 


rrTlf'^ + b = -/r'*/"’diVa;(T>e53). 


Proof: Insertion of (13.251) and (13.261) into ^3.271) yields 

= T>(£(V,n) - eS) + y^/^D{e{V^u) - eSi) 
+ yD{eiV^u) - eS2) - 


(5.28) 

(5.29) 


Using (I3.36p . we see from this equation and from p2.17p . (12.221) . (I3.29p that p5.28p holds. 
(I5.29P is an immediate consequence of (15.281) and (|2.16l) . (|2.21l) . (|3.28l) . ■ 


5.2 Asymptotic expansion of St + c(W 5 + — A.xS) 

In this section we compute the form of the residue 

{pX)^/‘^dtS + c(^ds\N{s{Vxu), S) + - y.^/^XAxS), (5.30) 

which is obtained when we either insert for (u, S) the inner expansion [u^^\ or the 

outer expansion of the asymptotic solution 

For functions {t, x) 1 —>■ w{t, x) defined in a neighborhood of T we write w{t, r], ^) = 
w{t,x) with X = rj + n(t,r])f^, as always. However, in the following computations this 
slight abuse of notation could lead to confusion when we consider derivatives with respect 
to t. To avoid this, we introduce the notations 

w\t{t,x) = w\t{t,ri,^) = drw{r,ri,^)\^^^, {dtw){t,r],^) = dtw{t,x). 

As introduced previously, for i = 0,1, 2 we write S''(t, p, C) = d^Siit, p, C) and S'/{t, rj, C) = 

d‘lSi{t,r],C). 
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Inner expansion We first compute ()5.30p for {u,S) = To this end we 

need 

Lemma 5.5 Let be the normal speed of the phase interface T{t) at p G L'{t), 

let be the operator defined in dMl), and let w be a function defined in a neighborhood 
ofT. Then we have 

dtw{t,x) = w\t{t,r],f) - ^{dtn){t,r]) ■ Vrjw{t,rj,^) - s^^^\t,r])d^w{t,T],f,). 

Proof: By definition, V^rc(t,ry,^) is a tangential vector to r(t). Lemma \2 ^ thus yields 

dtw{t,x) = dtw{t,ri,f) = + dth ■ V^w{t,r],^) + d^w{t,rj,f,)dtf 

= w\t{t,r],C) - C{dtn){t,r]) ■ \/rjw{t,rj,^) - s‘^^^\t,rj)d^w{t,r],^). 


This proves the lemma. ■ 

We apply this lemma to the function s[^'^ defined in (|3.50p to obtain 

dtS[^\t,x) = s['^^{t,'n,^) - ^{dtn){t,rj) ■ Vr,s['^\t,r],^) - s^f^\t,r])d^s['"\t,r],^). 

From this equation and from the asymptotic expansion (I2.44p of s^^^ we conclude for the 
first term in (15.300 that 

{p\f/^dtS[^\t,x) 

= (;.A)‘/28,(s„(_^) + {t ,,, + 8S2(t,,, i^)) 

= ~('So + S'l + P.S2) + pX^^'^Rst 

= —sqSq — p^^‘^{siSq + sq^i) + pRsti (5.31) 


with 


= {Si\t +p^/‘^S2\t} - fidtn) ■Vr,{Si +p^/‘^S2), 

Rs,ih,\t,g,^) = xy^Rg^-siS[-{so + p^/hi)S'2. (5.32) 

For the third term in ^5.301) we get from Taylor’s formula and from p3.50p 
^^'( 5 ;^)) = + ^"(5o)5i + /r'/2(^"(5o)52 + + pR^, (5.33) 

where 

Rf' = lf’'"{So){2SiS2 + p^^^S^) + ^^^^^\So + ^{p^^^Si+pS2)){Si + p^/^S2f. (5.34) 
with suitable 0 < < 1. Observe next that 

A^s[^\t,x) = d^s[^\t,p,C) - K{t,p,^)d^s[^\t,p,^) + Ar^s[^\t,p,^), (5.35) 
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where Ar^ = divr^Vr^ denotes the surface Laplacian and where K{t,r],^) is twice the 
mean curvature of the surface T^{t) at 77 G With the notation K'{t,r],0) = 

d^n{t,r],0) we obtain from Taylor’s formula 

K{t,ri,^) = K{t,ri,0) + d^K{t,ri,0)C+ ^d^K{t,ri,^*)f 

= nrit,r]) + (5.36) 

where = ^d^K{t,ri,^*) is the remainder term, with suitable between 0 and 

We insert (I3.50p and (|5.36l) into (I5.35P and obtain for the fourth term in (|5.3Up that 

+ S'{ + 77 ^/ 25 " 

- ^Acr + + fiXR^C^^ {Sq + - fiX^^^K{^)S2 + 

= + (-^i - A^/^Kr 5 (,) + 77^/2(52 - X^/\rS[ - Ak'C^q) + ( 5 . 37 ) 

where 

Ra = -Xk'CS[ - X^/\iOS'2 + AArJSi + ^^^^52) + X^^^R^C\S'o + (5.38) 

From (I5.3ip . (I5.16p . (15.331) and (j5.37p we obtain 

{^iXy/^dtS[^'> + c(55W(e(V,ui^)), S{^^) + 

= - '^o) + c(^"(5o)5i -S'I-e: [f]5o - <t(0) + {X^/\r - ^)S'o) 

+ 077^/2 ^V;"(5o)52 -S'^-e: [f]5i - <t(0) + {X^/\r - '-^)Si 

- X^^^(^a'{0)C + £ : De{a* <S)n + Vru*) 5o(7?)di?) 

+ 2'*^^^^('^o)<S'i + {Xk'C, + l^^St+c {...)) (5.39) 

where 

Rst+c(...) = Rst + c{-Rw + R^, - Ra), (5.40) 

with Rst, Rw, R^i and Ra given in p5.32p . p5.17p . ^5.341) and P5.38F respectively. 

Corollary 5.6 Let sq be given by ^2.4011 and assume that the functions Sq, Si and S 2 
satisfy the ordinary differential equations (I2.26p . (j2.27p . (I3.55p . with Fi, F 2 given by 
(I2.32p . (I3.56p . Assume moreover that the conditions (I2.28P - (I2.3ip and (13.601) . (13.611) 
hold. Then there is a constant K such that the interior expansion {u^i \ defined 

in p3.49p - p3.5ip satisfies for all {t, r],f) G Qinn ^ U and allO < ^ < po, 0 < A < Aq 

the inequality 

= [^^'\Rs,+c{-Rw + R^,-RA)\<K{^^^'‘'\\nyi\\ (5.41) 
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Proof: From (j2.40p we obtain 


C Cl 

and by (15.91) we have e : T^~'^ = (t( 0). After insertion of these two equations into ()5.39p . 
the latter equation takes the form 

= - S'^) + c(^"(5o)5i - S'{ - Fi) + - S'i - F^) 

+ iJ^Rst+c{...) = + c{-Rw + R^i - Ra)) ■ (5-42) 

Here we also used (I2.26p . (12.271) and (I3.55F Noting the inequalities (15.221) for ^ and C, the 
inequalities p4.3p - (14.51) . p4.11l) . (14.131) . (14.141) for Sq, Si, S 2 , and the inequality ^5.261) 
for the term Rt^, which appears in Ryy, we see by inspection of every term in (15.321) . 
(I5.17p . (I5.34P and (|5.38l) that the inequality 

l-^St + c{-Rw + R^, - Fa)I < K\ In^p (5.43) 

holds. To obtain inequality (|5.41l) we divide (I5.42p by and estimate the right 

hand side of the resulting equation using (|5.43l) . ■ 

Outer expansion Next we compute (I5.30p for {u,S) = (u 2 ^\ . Note first that 

pi.8p and p5.28p yield 

ds\N{eiV^ui^^), S!f^) = -e : (f + + ^f) + : DeSs . (5.44) 

Also, Taylor’s formula and (|3.26l) yield for a suitable 0 < '&{t, x) < 1 

^'(sj^)) = ^'(S) + + ^lS2 + 

+ ^^^'"(S)(//1/25i + fiS2 + 1x3/253)2 + + 1X52 + 1 x 3 / 253 )^ 

+ {S + ^/(ix^/'Fi + 11 F 2 + (ix'/25^ + ^§2 + 1 x 3 / 253 )^ 

= ix1/2<(5)Si + 11 (<(5)52 + lr’{S)Sf) 

+ 1x3/2(^"(5)53 + <'(5)<52 + (5)<3) + R., . ( 5 . 45 ) 

Here we used that '0'(5) = 0, by (|2.38p . Equations (I5.44p and p5.45p imply that in the 
domain Q \ T 

(ixA)1/29i 5< + c(asW(e(V.xx<), 5<) + - 1 x'/'AA,,5(^)) 

= c(-e:f+ <(<5i) 

+ cixi/2 ( - e : f + <(<52 + ^<'(«) 

/ 1 A ^/2 \ 

+ cix( - e : T + <(5)< + <'(5)« + - XAji + -^dtSij 

+ 1x3/^ Rst+c{...) ) (5.46) 
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where 


Rst+c{...) = c (e : DeS^ + R^, - XA^{S 2 + + X^^^dt{S 2 + (5.47) 

Here we used that the function S has the constant values 0 in 7 and 1 in 7 '. 

Corollary 5.7 Assume that the functions Si, S 2 and S 3 satisfy (|3.30p - (|3.32p . Then 
there is a constant K such that for all {t,x) ^ Q\T and all 0 < fi < fio, 0 < A < Aq 


dtS, 


iu) 


+ 


(/xA)V2 




< K 


Jf_ 

AV 2 ■ 

(5.48) 


Proof: By p3.30p ~ p3.32p . the brackets on the right hand side of equation p5.46p vanish. 
Therefore, if we divide the latter equation by (/rA)^/^, we obtain 


dtS^ 


in) 


+ 




^ 7? 

Rst+c(...)- 


AV 2 


5i and S 2 are given in (I3.36p . p3.37p . and the function S 3 is obtained by solving (I3.32p for 
this function. From these equations we see by our general regularity assumptions that 
ll(-5i,-52, 53)11 L°“(Q\r) < -^ 1 , with the constant Ki independent of fi. Using this, we see 
by inspection of every term in (|5.47l) that ||7?5j_|_c(...) ||Loo(Qyr) < with K independent 
of fi and A. This inequality and the equation above imply (I5.48h . ■ 


5.3 Auxiliary estimates needed in the matching region 


The following auxiliary estimates are needed to prove (12.511) and (12.5311 in the matching 

Q (/iA) 

match- 

Lemma 5.8 The functions , T!f'\ S^'^ defined in (|3.25D - P3.27D and u\^\ t[^\ 
s[^^ defined in (|3.49p - (|3.5ip satisfy 


II cin) _ qw II 

l|5l 

\D-{S[^'> - si^'>)\\ ,,,, 

^ '^^match 


:(m) I 


|q| 3-|c 

2 /i 2 


3t{s[^'^ - S^2^^] 

iu) ,,(m)i 




11^1 

^ '^^match-' 

IIT^Cm) rpiu) II 

II 1 2 


< A//3/2|ln^|2 

< KX 

< KX-R^fi, 

< KXR^yi^/^\lnyi\, 

< Kn, 

< Kn, 


1 < |a| < 2, 


(5.49) 

(5.50) 

(5.51) 

(5.52) 

(5.53) 

(5.54) 


for all fj, G (0,/io]; A € (0,Ao]. Here a denotes a multi-index and K denotes a positive 
constant, which does not necessarily have the same value in the six estimates. 

Proof: Since the proofs of these estimates are long and technical, we present here only 
the proofs of ^5.4911 and p5.52p . The proofs of the estimates p5.50p , p5.5ip and p5.53p run 
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along the same lines. (15.541) is an immediate consequence of the definitions (j.S.27p . (j8.5ip 
of and , and of the estimates (15.531) , (I5.49P . 

In this proof we mostly drop the arguments t and r] to simplify the notation. As usual 
we write C = need that for G Q^atch inequalities 


In^l 


< 


(^A)V2 


= IC|<3 


lnn\ 


(5.55) 


hold, by definition of Q^atch ™ (!2.34p . 

We begin with the proof of (I5.49P . By dehnition of and in (I3.50p and (13.261) 
we have 


< \So-S\ + fi^/^\Si+n^/^S2-Si-iJ.^/^S2\ + iJ.^^^\S3\. (5.56) 

To estimate the first term on the right hand side note that since S{t,x) = S{^) = l'''(0) 
relations (IO]l . (lOll . and (|5.55l) imply for {t,r],^) € Qj^^'ateh 






(5.57) 


To estimate the second term on the right hand side of p5.56p we introduce the notations 


Pi(C) 


P2{0 


ip"( 0 ) ’ 
£:T’(+) 

b"(i) ’ 
1 


1 

b"(i) 


C<o, 





C>o, 





: - 

W'(0) 1 

1 b"(o )) 

+ As (t'(O)C^ , 

C<o, 

: - 

b"'(i) 1 

(e-.fW\ 
K W(i) ) 

+ A 2 (t^( 0 )^ 



“h A 2 £ 

: De{a* 

®n + Vrn*)c), 

C>o. 


(5.58) 


(5.59) 


By definition of the functions pi, p 2 in (I4.10|) and (I3.59|) . and by definition of ip, p± in 
P3.57p and ^3.581) . we have 


Pi{C) - Pi{C) = (1 - V^(-C) - <7>(C))ft(C) = 0, for i = 1, 2 and |C| > 2. 

We can therefore choose a suitable constant K such that \pi{C) — PiiC)\ < holds 

for i = 1,2 and all C € M. From this inequality and from the estimates (|4.12l) . (14.151) we 

conclude that 

\SiiC) -Pi\< |5i(C) - Pi(01 + |pi(0 - Pil < (i^2 + i^)e-“l^l, (5.60) 

1^2(0 - P2| < 1^2(0 - P2(0I + |P2(0 - P2| < (i^5(l + ICl) + (5.61) 

for C G M. 

Now we proceed to estimate the second term on the right hand side of (|5.56l) . We 
insert the functions pi and p 2 into this term, use the expressions for Si and S 2 given in 
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< 


Pi + P^^‘^P2 - f 


(j3.36p . (j3.37p . and employ the triangle inequality to obtain 

\Si+p^/^S2-Si-p^^^S2\ 

' e:f j)"'{S) ^ e : f 
ip"{S) ' 2V>"(5) ^ 

+ 1^1 - pi\ + /iV2|52 _ P2I = |/,| + I/2I + I/3I. (5.62) 

By (I5.60[l . (15.Gill and (|5.55p we have for € Q^atch 

1^21 + 141 = lSi(t,p,C) - Pi(i,'r/,C)l + P^^^lS 2 (t,p,C) - P 2 (t,p,C)\ 

< (K 2 + + p^^^K5(1 + Id) + K)e-“lf 

< (^C7i+^^/2(272(^i + 3H^^^e-f|inH <C3/i3/2. (5,63) 

To find an estimate for |/i| note that the definitions of pi, p 2 in (I5.58p . (I5.59P yield for 
(t, p, 0 € with C > 0 


h = Pi 


(^A)V2 


+ P'^P2 




(^A)V2/ d"(l) 


e:r(0 i/e:T(0 ii'" {!) fe : T{0\^ 

- p2 


d"(l) 2d"(l) ^ V’"(l) 


d"(l) 


: f W - e : f(0 + (T'(O)e + e : De{a* ®n + Vp^*)? 


(5.64) 


d"(l) 


(j5.2l) and (I5.55P together yield 

|e : fW + cj'(0)e + e : De{a* ®n + Vtu)^ - e : f (e)| 

= < C'4MA|ln/ip. (5.65) 

Since = T(0+), = T'(0+), the mean value theorem and p5.55p imply 


P 


1/2 




= p^/^\R{^)^\ < Cs/iA^/^lln/rl, (5.66) 
where the remainder term R belongs to L°^{Us). Combination of (I5.64p - p5.66p yields 
for (t, p, 0 € Qldateh with ^ > 0 that 

|.fi| < p{C 4 ,\\ In/ip + CsA^/^l ln/i|) < Cq A^^^j In/ip/i. (5.67) 

From the definitions of pi, p 2 in p5.58p and p5.59p we see that for (t,p,^) G Q^atch with 
^ < 0 the term R takes the form 

1 


d"(o) 


Ii = -.R— : f (-) - e : f{() + a'{0)( 

+ /iV2(e:r(-)-e:f(0)-p 


1/2 


e : f(-)\2 

d"(o) 


^ ^ d^'(o) 
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Using (15.3p instead of (j5.2|) . we see as above that the estimate (j5.67p also holds in this 
case, whence the estimate (15.671) is valid for all 0 ^ ^matih- 

To finish the proof of (j5.49j) . we combine (I5.56h with (I5.57h . (j5.62j) . (I5.63h and (j5.67j) 
to obtain the estimate 

which implies (I5.49h . 

Next we prove (j5.52p . Prom (I3.25|) and (13.381) . (j3.39p we conclude for {t,x) € Us that 


{t, x) = u{t, x) + x) + x) 




e : f (+) , e : f (") 

-? “I--? 

V’"(1) V^"(0) 


+ v(t, x) + x) + fj,u{t, x). 

Combination of this equation with (I3.49h and insertion of (I3.52h - (I3.54h yields with 


c = 


« 




that 


uf (t,x) — U 2 {t,x) 






= So W d7? - C+) + /uA 

J — 00 ' 

e : f (+) 


e : e : 


-C + 

C . i-d 


i ) 


+ /rA^'^^u’ 


— 00 

C 


i 


e ■ T^~'^ 

l+(r?)-;- 1 


= (/iA)i/2ji(c)+ ^i{J2ic) + J3ic) - m)- 


a 17 

So(r9i)d??i -r?+) di} 

-00 ' 

did — fJ,u{C) 


(5.68) 


To estimate the right hand side we use the boundary condition (I3.35P for u, note that 
by (I5.58P the equation pi{C) = 7 ^^ holds for C > 0, and employ the inequalities (|5.60p . 

(|4.8n to compute for {t,r],^) G ^ > 0 

/ OO /*l 9 

(/ So{idi)did,-^+)d^ 

-OO J — 00 ' 

+j3(c)-Avvr--'^ 

Jo 


ril) 


did 


roo , [•'& \ foo 

Xa* (y So(i?i)^9i-^9+)di? + A^/Vy^ (Si(i?) - pi(C)) 

roo poo 

<A|a*|/ —e-“’’dr9 + A^/Vl / {K 2 + d-d 

Jc « Jc 


dd 


,Ki 


{X\a*\^ + ^ ||r(*||ioo(r))e-<. (5.69) 
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The mean value theorem implies for ^ > 0 


u{t, x) = u{t, r], C) = u{t, T], 0+) + Ru{t, r/, 0?- 

Since ■u(t,r/, 0+) = we infer from this equation and from ()5.69p . (15.551) for all 

(t, rj, 0 G Qmatch '^ith ^ > 0 that 

1^2(0 + JsiO - m\ < \M0 + MO - + |n(0 - u(+)| 

< C'iA^/^(||o*||ioo(r) + ||tt*||ioo(r))e + \\Ru\\l°°{Us)^ 

< C'iA^/^(||a*||ioo(r) + ||tt*||ioo(r))e + ||i?ii||£,oo(^^)-(//A)^/^| ln^| 

< C' 2 A^/^/x^/^| ln/x|. (5.70) 


Using the boundary condition (I3.34p for u instead of (I3.35p . we see by the analogous 
computation that (I5.7nh also holds for {t,r],^) G Q^atch ^ < 0. 

Now use p5.70p to estimate the second term on the right hand side of p5.68p . The first 
term is estimated by p4.8p . Because p5.55p yields we obtain for 


{t,v,0 e Q 


ifiX) 

match 


\O0\t,x) - 02^\t,x)\ < f SoiOd'd-C^ +/^|-^2(C) + <^3(0 - «(0| 

J —oo 

< (/rA)^/^(max |u*|) —In/r| < In/i|, 


which implies (I5.52h . 


5.4 End of the proof of Theorem 12.31 

To complete the proof of Theorem 12.31 note first that p3.22p , p3.23p imply 


0-^^ I = vM 
lao ^ 


ao’ 


o g{R) I _ o e(M) , 
^rign‘^ ~ '-'nan‘^2 


lao’ 


cIm), _ C'(/^)| 

Iq(mA) 1 Iq(('^)' 

^ inn ^ inn 


Therefore p2.56p follows from the definition of in p3.50p . equation p2.45p follows from 
p2.28p . p2.3ip . p3.60p . and p2.49p is a consequence of the definition of ™ ^3.251) 

and of P2.20I) . (|2.25l) . (|3.33l) . Moreover, the estimate (|2.55l) for the right hand side 
of (|2.50l) follows from the definition of in (|3.26l) and from dnan<S| = 0. This 

last equation holds, since by assumption r(t) C fl, which implies that S{t) is identically 
equal to 0 or 1 in a neighborhood of dfl. 

(12.471) follows from the definition of in (|3.24l) : equation (|2.57l) is an immediate 
consequence of (|4.14l) . 

It remains to verify the estimates p2.5ip - p2.54p for the right hand sides 
of the equations p2.46l) and p2.48l) . To this end we put together all the estimates derived 
in Sections 15.11 - 15.31 We start with the proof of (I2.5ip and (|2.52l) . 

Equation (I3.22p yields 
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We insert this equation into (I3.24|) and use (I3.23|) and (13.271) . (I3.5ip to obtain 


+ T!f\l - cl>^x) + De{{u[^^ - u[^^) ® (5.71) 

The function (f)^\ defined in (|3.2ip is independent of (t,T/). The decomposition (13. 5p of 
the gradient thus yields 

with the unit normal vector n = n{t,r]) to r^(t) at r/ G r^(t). We write 

^ ^ ^3(/iA)V2|ln^|^’ ^ ^3(pA)V2|ln/ip’ 

by a slight abuse of notation. With p5.7ip and (15.721) we compute 
diVa^r^'^) + b = (diVa^r^^^^ + b)^^A + (diVa;r2^'^^ + b)(l - ())^a) 

+ ((r*"* -Tj'‘')n + div,Dg((iii'‘' 

+ (0e((„<'‘> - u'"’) ® „))„ (5,73) 

Inequality (I2.52p is an immediate consequence of this equation and of (15.291) , since (/>^a = 0 
in (5out\ and p2.5ip is obtained by estimating the right hand side of (15.731) using the 
obvious inequality 

|diva:lle((u^^^ — ^ 2 ^^) (8> n)| < C[\Vx{u[^^ — U 2 ^^)| + \u[^^ — U 2 ^^\) (5.74) 

and the equation and inequalities ^5.211) . p5.29l) . ^5.521) - p5.54p . 

We next proof (I2.53P and (I2.54p . The inequality (|2.54l) follows immediately from 
(I5.48p . since = 0 on Qout^ which by (|3.22l) and (|3.23p implies ) = 

{4^\si^^)on Qi^4'>CQ\r. 

It remains to verify (I2.53p . Since \Ns{e, S) = —e : D{e — eS), by pi.Sp . it follows from 
(IMP . (I33TP . and (IMP that 

\Ns{e{VxU^^^), 5(^)) - \Ns{e{Vxv4^), s[^^) = -e : (T(^) - 

= -£ : (t4^ - t4^)(1 - </>ma) - e : De((J4 “ 4""^) ® 

The mean value theorem and (13.231) imply 
_ ^'(5^) = 

for a suitable 0 < i!}{t,x) < 1, and 

AxS^f^'> - AxS[^'> = AxiS^'^ - - ^^x) + - s4'>) • V,</.^a 

+ - s4^)Axct>^x- (5.75) 
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The right hand sides of the last three equations vanish on the set since (/>^a = 1 

on this set. On the set right hand sides of these equations can be estimated 

using (15.491) . ()5.50l) . ()5.52p . (15.541) . ()5.72p . In the estimation of (|5.75p we also note that 
since (/>^a is independent of analogous to (I5.35p the equation 






3(A/x)i/2| ln|u|*^'^'^ 


( 2 a n 2 

V3(A/r)^/2| ln/r| / 


holds, with twice the mean curvature K{t,r],^) of the surface r^(t) at rj £ r^(t). Together 
we obtain that on ^ U Q^^tch inequality 


< K[^ + In ^1^ + < K^\ In^up (5.76) 


holds. Similarly, (13.231) implies 

- dts[^^ = - ct>^x) + 

The right hand side of this equation vanishes on To estimate the right hand side 

on the set inequalities p5.49p . p5.5ip and the equation 




2asi^) 


^4^fj,xi 


3(/iA)^/^| In fj,\ 

which follows from ^3.211) and Lemma 15.51 The result is 

< iLA-VVl ln^^l on U 

By combination of ^5.411) . p5.76p and (15.771) we see that the inequality 


(5.77) 


< 


qM 

X^l 


+ K\ ln/i| + 


(^A)V 2 


Kijl\ In^l" 


< K 


© ' + ^ ' |ln/i|2 

holds on the set U Q^atch- proves (12.531) and completes the proof of Theo¬ 


rem 
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6 Proof of Theorem 12.81 


This section contains the proof of Theorem [221 (u, T, S) denotes the asymptotic solution 
constructed in Theorem 12.31 and (uacj Tacj 5'ac) denotes the exact solution of (jl.ip - 
(11.31) . ()2.67l) - ()2.69l) . For the proof we need a lemma and a theorem, which we state 
first. 


Lemma 6.1 For all x G T{t) the propagation speeds sac <^nd s satisfy 
SAcii,x) - s{i,x) = |y ^ • i'^Aciix) -T{i,x))y (6.1) 

where is the right hand side of equation (I2.48P . 


Proof: Since the manifold T is a level set of S and since by (12.701) the manifold Tag 
is a level set of SaCj h follows that {dtS{t,x),VxS{t,x)) and ((9t5Ac(L 3;), V 3 ; 5 ac(L 3;)) 
are normal vectors to the respective manifolds at {t,x). Moreover, (I2.69P implies that 
Vx^acCL 3:) = VxS{i,x). From (j2.9p we thus infer that 


s{t, x) = 


SAc(i, x) 


—dtS^i,x) —dtS{i,x) 

VxS{i,x) ■ n{i,x) |Va,5(t,x)|’ 
_ -dtSAcji, x) _ -dtSAc(lx) 
~ |V,5ac(T 3:)| “ |V,5(Lx)| ■ 


( 6 . 2 ) 

(6.3) 


For brevity we do not write the argument (t, x) in the following computation. In (|6.2I) 
we eliminate dtS with the help of (|2.48p . and in (j6.3p we replace dtS ac by the right hand 
side of (jl.3p . Together with another application of (j2.69p this results in 


«AC - S = ^^^3^y^^^^((95W(e(V.UAc),5Ac) + 

- {dsVJ[e{Vxu), S) + - /ri/2AA,5)) + 

= ^ “ (/xA)V2|v,5| • ^AC - e : t). 

which is (j6.ip . In the last step we used that by (11.81) and (12.471) we have 

95W(e(V,nAc), Sac) = -e ■ Tag, and ds\N{e{Vxu), S) =-e : T. 

The proof of Lemma l6.II is complete. ■ 

Theorem 6.2 Suppose that the order of differentiability of ip, F, u, u, b, is higher by two 
than required in Theorem \2.^ Assume that the principal eurvatures of the 

regular -manifold T{i) = F(^'^)(t) are bounded, uniformly with respect to p, G (0,/tq] 
and A G (0,Ao], and that there is an open subset Q' CC P and <5 > 0 such that the 
neighborhood o/r(^'^)(£) defined in (12.31) satisfies C Q'. Then there is a 

constant such that for all p G (0, po] and all A G (0, Aq] 

||Tac - ?'llL2(r(t)) ^ In^lV- (6.4) 
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We postpone the proof of this theorem and hrst finish the proof of Theorem 12.81 


End of the proof of Theorem 12.81 Bv (12.7p and (12.561) we have for x = (r/, 0) € T{i) 
that 

VxS{i, x) = n{dnS){i, x) + VrS{i, x) 

+ r/, 0) + ^Vr52(t, r], 0). (6.5) 


(14.11) implies 

5(,(0) = ^2^(So(0)) = ^2^(l/2) > 0, 
whence, from (j6.5p for /r € (0, ^o] and A € (0 , Aq] with fiQ sufficiently small, 

\VxS{i,x)\ > - h|5c5’2(t, r?, 0)|) 

- H^/'^\VrSi{i,r],0)\ - n\VrS2{i,r],0)\ > ( 6 . 6 ) 

Combination of (16.ip with the inequalities p2.53p . P6.4I) and (|6.6I) yields 


\\sAcii) - s(i0llL2(r(£)) 

1 


< 


minr(£) iV^.S'Pt)! 

2 (/rA)V 2 


2^" ^(^lllL2(r(f)) + (^a)V 2 II^Ac(i) - T{i)\\^2(Y'(i))^ 
(|ln/^P(0 ' i^3meas(r(t))^/2 ^-^-^I^Kglln/rlV) 

< Kq\ In/rp/r + Kj\ ln^|^/r. 


< 


2 V'(l/ 2 ) 


p2.75p follows from this estimate. The proof of Theorem 12.81 is complete. ■ 

Proof of Theorem 16.2t Note that the function (uac(^)) ^Ac(c) solves the equations 
pi.ip . pi.2p in O with 5 ac(^) = S{i), by the initial condition p2.69p . Moreover, p2.67l) 
holds for UAc(c- From the equations (|2.46l) . (12.471) . (12.491) we thus conclude that the 
difference (uac ~ u, Tag — T) satisfies 

-div,(rAc-r)(tl = -/J^^^(t), (6.7) 

{TAc-T){i) = De{Vx{uAC - u)ii)), (6.8) 

{uAC —'u){i,x) = 0, X G (9n. (6.9) 


This is the Dirichlet boundary value problem for the elliptic system of elasticity theory in 
the domain Q. It suggests itself to derive the inequality (16.4p by using the L^-regularity 
theory of elliptic systems, which allows to estimate the norm UTac ~ F||^ 2 (r(f;)) by the 

L^-norm of the right hand side of p6.7p . To apply this theory directly we would 

need that the L^-norm of decays to zero for ^ ^ 0 uniformly with respect to A. 
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However, the relation meas(Q|^^^(t) U < Ci(/iA)^/^| ln/x|, which follows from 

(j2.34p . and the estimates ()2.5ip . p2.52p yield 


3/4 


The right hand side does not decay to zero for ^ ^ oo uniformly with respect to A, 
but blows up for A —?■ 0. Therefore direct application of the L^-regularity theory is 
not possible. Before giving the detailed proof of (|6.4p we sketch how to circumvent this 
difficulty. 

Set 

A(//) =lnp|, (6.10) 

a 

where a > 0 is the constant dehned in (j2.38p . By (I2.34p we have 

Qi™’(«) u QL'llhffi = {(x.O 6 I l«l < . (6.11) 


Define 




I), i-€(3LT(«)uq|;‘1A(<) 

6 . 


( 6 . 12 ) 


hence For x = x{i,r],^) G Us(i) we write as usual f[i^\x) = 

From (|2.5ip and (|6.11l) we obtain for rj € F(t) and —A(/r)A^/^ < ? < A(/i)A^/^ 

that 

l/iT’(il.0l<|lni<t(9‘'bl'i. (6.13) 

Define : F (t) ^ M by 

iv) = Ai/2/Jf ^ {ri, dQ. (6.14) 

(I6.13P and (16.101) together imply that 

< 2A(/r)A^/^| In^l^ ^ Ki = In/rp, (6.15) 

for all r] € r(t). Examination of the boundary value problem (16.71) - (16.91) suggests 
that for /i hxed and A —>■ 0 the solution (uac “ u, Tag — T) converges to the solution 
(u*, r*) : D ^ X of the transmission problem 


div3;T* 

= 0, 


(6.16) 

n 

= De{V, 

cti* ) ? 

(6.17) 

[T*]n 

- 6^^'> 

on r(t), 

(6.18) 

[ti*] 

= 0, 

on r(t), 

(6.19) 

tt*(x) 

= 0, 

X G 5D, 

(6.20) 
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where 6i^\r]) = si^^\r]) for 77 G r(t). If this limit exists, it follows from (j6.15ll 

that 

a 

This implies that the solution (u*,T*) will be bounded by C/r|ln|Up with a suitable 
constant C, and this limit behavior suggests that though the L^-norm of blows 

up for A ^ 0, the solution (uac “ u,Tac — T){i) of (I6.7p - (I6.9p is bounded by C^\ lnfi\^ 
with C independent of A. The reason for the blow up of ||/i^^^(t)||L 2 (r 2 ) for A ^ 0 is 
therefore not that the norm of the solution (uac ~ u, Tag — T) (t) would blow up, but that 
the solution looses regularity in a neighborhood of the surface T(t), which is shown by 
the equation (|6.18l) for the limit solution. This equation implies that T* does not belong 
to the Sobolev space 

In the following proof we do not study the limit (u*,T*). Instead, based on the idea 
of the behavior of the regularity of {uac — u, Tag — T){t), we decompose this function in 
the form 

{uac - u,Tag - T){i) = 

where = is bounded by C/illn/rp, uniformly with respect to A, 

and for A ^ 0 has the same regularity behavior as {uac — u,Tac — T), but otherwise 
does not approximate {uac — u,Tac — T). The construction is such that the difference 
ri^^) = ri^^) = (uac — u, Tac — T) — {u^^\T^^'>) does not loose its regularity 

for A ^ 0. Hence, we can use the standard L^-theory for elliptic equations to show that 
also {u^^\t^^^) is bounded by C^\ In^up independently of A. 

To construct {u^^\T^^^) let Us{i) be the neighborhood of r(t) defined in (12..Ih and 
let € C'^((—(5, (5)) be a function satisfying 

(t>,{i) = 1, -8/2 <i<8/2. (6.21) 

We set 



X = x{i,r],C) G Us{i), 

(6.22) 

lo, 

X G Q\Us{i), 


tW(x) = De{S/,u^^\x)), 

X G XI, 

(6.23) 


where the function Q i-)- V{X,r]X) ■ [— Xl/2 ) ;^l/2] is constructed as follows: We use 

the notations V = d(Y, V = d‘^V. In the interval [—A( 7 r), A(/i)] the function V is the 
solution of the boundary value problem 

(^T)e(H"(A,r/,C) <8 )n))n = A^/Vii^^(^, -A(7r) < C < A(Ai), (6.24) 

V{A,rj,±A{^)) = 0, (6.25) 

where n = n{i]} is the unit normal vector to r(t) at rj G r(t). The equation (I6.24p is a 
second order linear system of ordinary differential equations for the three components of 
V, which can be written in the form 

BV" = (6.26) 
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with a 3 X 3-matrix B = B{r]) defined by the equation 

Buj = [De{u} 0 n))n, (6.27) 

which must hold for all a; E The matrix B is symmetric and positive definite 

uniformly with respect to r]. To see this, note that since the elasticity tensor D : ^ 

is a linear, symmetric, positive definite mapping, we compute for ui,(jJ 2 E 


(Bui) • 002 = y{Ds{uJi (?) n))nj • L02 

= {i 02 n) : De{iOi n) = e{i 02 < 2 ) n) : Ds{loi 0 n) 

= [De{uj2 ( 8 ) n)) : e{uji (g) n) = (^[De{uj2 <S> n))n^ ■ coi = {BUJ2) ■ wi- 

This shows that B is symmetric. For w E we have with a suitable constant Cq > 0, 
which only depends on D but is independent of r], that 

(Boo) ■ 00 = e{oo (g) n) : De{oo 0 n) > C'o|e(w 0 n)f > -^\oo\^, 

hence B is positive definite uniformly with respect to r] E T(t). 

Therefore the boundary value problem (16.241) . (16.251) has a unique solution V on 
[—A(//), A(/i)]. To extend C to all of [—we continue V to the 

intervals ( — — A(/r)) and (A(/r), by affine functions: 


v{x,vX) 


(c + A(m))F'(A, V, -A(^)), -3^ < C < -A(m), 
(c - A(/r)) F' (A, r/, A(/x)), A(/x) < C < • 


By this extension, ( 1 -^ V{X,t],() is continuously differentiable at C = iA(/r). For 
X = x{t,T],^) E Us we use the notation 

F(A,x) = F(A,r?,^). 

In the remaining part of the proof of Theorem 16.21 we need the following lemma, which 
we prove first. 


Lemma 6.3 There are constants C'i,...,C '4 such that for all E (0,/ro], A E (0, Aq], 
(?7, C) G r(t) X ( — CL'iT'd X E Us{i) the estimates 


\Vi^V'{X,v,C)\ < Ci|ln/r|V, j = 0 ,l, 2 , (6.29) 

|Ai/2v);F(A,r/,C)| < C2|ln/r|V, j = 0,l,2, (6.30) 

|A^/2v^F(A,x)| < CsIln^lV, (6.31) 

|A^/25,,Vr,F(A,x)| < C4|ln^|V, A: = l,...,3, (6.32) 

hold. Moreover, there is a function —>■ and a constant C 5 such that 

defined in (|6.23p satisfies 

div.rW = + (6.33) 

with 

|5^^^)(x)|<C5|ln/r|V (6.34) 


for all X & Q and all fi E {0,fio], X E (0,Ao]. 
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Proof: In the following computations we drop the arguments A and rj. Integration of 
(j6.26l) yields 


BV'{C)= +W'(-A(^)), (6.35) 

J 

BV{C)= r AV2/}f^(AV2^?)dMi?i + (C + A(/r))W'(-A(/r)), (6.36) 

J —A(^) J —A{p) 

where we used the boundary condition (j6.25p to get the second equation. Since V (A(^)) = 
0, the relations (|6.36l) and (|6.13p together yield 


2A(/r)|W'(-A(/x))| 


hence, by (I6.10p . 



r'&i 




-A(/i) 




/■i?i 


< / 

/ 


J — A(^) o 

'-a(a) 



2A(/r)^|ln/i|^//^/^A:i , 


|W'(-A(m))| < A(M)|ln^|V/'^i = -ATilln^lV 

a 

Since B = B{r]) is positive definite uniformly with respect to rj, this inequality implies 
the estimate (I6.29P for j = 0 and — A(/r) < C < A(/r). Since by definition ()6.28l] we have 
^^(C) = ^^(~A(/r)) for C < —A(;u) and V'{C) = V'{A{p) for A(/i) < (, the estimate 
(j6.29l) with j = 0 holds also for the values of ( outside of the interval [—A(;u), A(/r)]. 

To prove (j6.3np for j = 0 we use that V{—A{fi)) = 0. By integration we thus obtain 
from (j6.29p for C € y^] that 


inoi = 




V'{^)d^ < 1C + A(;u)|Ci| In^lV < (A-^/2j + A(/i))Ci| In^ujV, 


which implies (|6.3UI) for j = 0. 

To verify ()6.29l) and (j6.30p for j = 1,2 we differentiate the differential equation (j6.26l) 
and the boundary condition (I6.25P with respect to t]. For j = 1 we obtain the differential 
equation 

B{v){dr„vr = - 5,,B(7?)B(7?)-Vjf ^), 

and a similar equation for j = 2. We then use the estimate 

|V(;/Jf^(^?,C)l = |V^,/}'^")(r?,C)| < |ln^|2(^)'^V, j = 1,2. 

This estimate is obtained by differentiation with respect to rj of the asymptotic expansions 
in Section 15 .1 1 leading to Corollary 15.31 Under the regularity assumptions in Theorem 16.21 
these derivatives exist. With this estimate we can employ the same arguments as above 
for the case j = 0 to derive (|6.29l) and (|6.30p for j = 1,2. 

To prove (16.311) we use the decomposition (13.611 of the gradient and (13.131) to compute 


V,U(A, x) = d^V{X, V, ® n + Vr^U(A, rj, ^) = A-^U' > 


n + (V^U)A(£,r/,C). 
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The right hand side is estimated by (16.291) and (I6.3n|) to obtain (j6..Sip . The estimate 
(j6.32l) is obtained from ()6.29p and (j6.30l) by similar decompositions. 

To prove (j6.33p . (j6.34p note that by (I6.22p . (I6.23p we have = De{\/xU^^^) = 
De{Vxi\^^'^V4>^)). Using p.6p and (|3.7p we therefore obtain by a similar computation 
as in (j5.73p that 

div,r(^) = diYxDe(Vx{X^^^V{X,7],^)M0)) 

= [De{V'' 0 n))n + divr^iAe(U' 0 n)'j (/>* 

+ A^/2 (^Be(d^Vr^V))n + divr^7?e(Vr^U))</>* 

+ ^(De(X^'^^VxV))n + diVxL>£(X^^^V 0n)'j(/>'^ 

+ ^T)e(A^/^U 0 n)^n4>'l 

= (6.37) 

In the last step we used the differential equation (I6.24h and noted that for ^ G 
([—<5, — A(//)A^/^] U [A(^)A^/^, d]) we have V”[X,r], = 0, by definition of V for such 

values of ^ in (16.281) . We also used that 4’*{0 = 1 for ^ G [—A(/i)A^/^, A(|u)A^/^], which 
follows from p6.2ip and ph.lOp . since we have chosen fiQ and Aq small enough such that 
A(/r)Ai/2 < 5/2 for all 0 < /r < Ho and 0 < A < Aq- 

The function is the sum of terms number 2 to 7 in the middle expression of 

equation p6.37l) . If we examine everyone of these six terms and apply (|6.29l) - (|6.32l) and 
also note that the functions (/>*, and cf)" are bounded independently of fi and A and 
vanish outside of Us(t), which follows from (/>* G C^{{—5, 5)), we see that p6.34p holds 
for . This completes the proof of Lemma 16.31 ■ 

To conclude the proof of Theorem 16.21 let be the solution of the boundary 

value problem 

-div^r?) = (6.38) 

= Ds{Vxui^'>), (6.39) 

ui^\x) = 0, X G dQ. (6.40) 

From these equations and from (16.221) . p6.23p . p6.33p we see that the function + 
) satisfies 

-div,(T(^) + ri^)) = 

(rw + ri^)) = D£(v^(uW + nl^))), 

+ ul^^)(x) = 0, X G 

hence is equal to the unique solution of the boundary value 

problem p6.7p - p6.9p . which means that = {uac — u,Tac — u){i)- 

Consequently, we have 

\\Tac - TWmm) < r^"^llL^(r(f)) + ri"^llL^(r(f))- (6-41) 
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To estimate |lL 2 (r(t)) we can use the theory of interior regularity for the elliptic 
boundary value problem (|6.38p - (|6.40p . By this theory there is a constant C such that 

where is the subdomain of 14 introduced in 
Theorem 16.21 hence by the Sobolev embedding theorem and by (I6.39P , 

||ri")|li2(r(,~)) < (6.42) 

where by our assumptions on in Theorem 16.21 the constants Ci,C 2 can be cho¬ 
sen independently of g and A. By definition of in (|6.12l) and by (12.521) we have 

|/^ 2 '^^(x)| < for all x € Q. From this inequality, from (|6.34p and from (16.421) we 

conclude that 

<C2{C5\]niifii + dxj (6.43) 

From (I6.23p . (16.221) and from the inequalities p6.30p . (I6.3ip we infer that 

|rW(x)| < c|v,uW(x)| = c\v4x^/^v{x,x)M0)\ 

= C\{X^/^V2,ViX,x))M0 + A^/V(A,x) ® {ncj)',{0)\ < 

whence 

Combination of this inequality with p6.41l) and (I6.43P yields p6.4p . The proof of Theo¬ 
rem [62] is complete. ■ 
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